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Abstract 

We consider the solution to a stochastic differential equation with a drift function which depends 
smoothly on some real parameter A, and admitting a unique invariant measure for any value of A 
around A = 0. Our aim is to compute the derivative with respect to A of averages with respect 
to the invariant measure, at A = 0. We analyze a numerical method which consists in simulating 
the process at A = 0 together with its derivative with respect to A on long time horizon. We give 
sufficient conditions implying uniform-in-time square integrability of this derivative. This allows 
in particular to compute efficiently the derivative with respect to A of the mean of an observable 
through Monte Carlo simulations. 

Keywords: Stochastic differential equations, invariant measure, variance reduction, Feynman-Kac 
formulae. 


1 Introduction 

We are interested in methods to compute the response of a Brownian dynamics to an infinitesimal 
change of a parameter A € R. More precisely, we consider the dynamics in R'*: 

r dA7 = Fx{X,^)dt + V2dWt, 

I X^=Xo, ^ 

for A € M close to 0, where {Wt)t>o is a standard d-dimensional Brownian motion independent of 
Xo ^ Note that neither the initial condition Xq nor the Brownian motion depend on A. The 
family of vector fields F\ ^ is indexed by a real parameter A. We assume that when A = 0, 
the vector field derives from some potential energy 1/ : ^ M, namely 


Fo = -VV, 
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where V denotes the gradient operator with respect to the space variables. For A close to zero, one 
can think of (Xt')t>o as a physical system undergoing a potential energy V to which one applies 
an external force Xd^Fx. Here and in all the following, the notation 9° denotes the derivative with 
respect to A computed at A = 0. 

Concerning the potential V, we assume that the following assumption holds. 

Assumption (Pot). The function V satisfies the following assumptions: 

(i) y : ^ R is a function such that x i—>■ V^V{x) is locally Lipschitz. 

(a) [ e~^^^^dx = l and f |Vy|^(a;)e~'^^^^da; < oo. 

jRd. Jgd 

(Hi) Pathwise existence and uniqueness hold for the process {Xt)t>o- 

Since W is assumed to be locally Lipschitz, pathwise uniqueness is automatically ensured. 
Pathwise existence is ensured for instance as soon as there exists a finite constant C such that for 
all X G R'*, XVix) ■x<C{l + |a;p). 

At A = 0, the dynamics © is of the following gradient form 

AX°t = -VV{X?)dt + V2dWt. (2) 

Under the above assumptions, it can be checked that e~^^^^dx is the unique invariant probability 
measure (see Lemma [18] below) denoted in the following: 

dvro = 


Then, from classical results in ergodic theory, for any / in L^(7ro), almost surely. 


\im J r f{X°)ds = f fdTTo. (3) 

Let us now introduce the assumptions we need on the drift (FA)AeiR. 

Assumption (Drift). There exists Ao > 0 such that, for all A G [0, Ao], 

(i) The function F\ — Fo is bounded by C\ for some constant C not depending on x. Moreover, 

as A —>■ 0"^, converges locally uniformly for a; G R*^ to some limit dxFx. Note 

that dxFx is bounded by C. 

(ii) The function x i—>■ Fx{x) is locally Lipschitz. The function x dxFx{x) is differentiable on R'* 
and V • dxFx is in L^(7ro). 

Under these assumptions, we will show (see Lemma [18] below) that the dynamics ITJ is ergodic 
with respect to a probability measure tvx: for any / G L^(7rA), almost surely, 

lim 7 [ f{Xs)ds= [ fdnx. (4) 

t->-00 t Jq J^d 

The aim of this paper is to study the quantity: for a given observable /, 


/gd /d^TA Jgd fd^o 
A —^0 A 



(5) 


In particular, we will exhibit sufficient conditions for the existence of this derivative, derive various 
explicit formulae for this quantity, and discuss numerical techniques in order to approximate it. 

The estimation of derivatives of the form ([SJ is useful in various applications, in particular in 
molecular simulations (see for example the recent work [34]): optimization procedure to fit a force 
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field to some observations, study of phase transitions, estimate of forces in Variational Monte Carlo 
methods (see mi or computation of transport coefficients. Transport coefficients are computed as 
the ratio of the magnitude of the response of the system submitted to a perturbation in its steady- 
state to the magnitude of the perturbation. These coefficients are related to macroscopic properties 
of the system through fluctuation dissipation theorems [nun]- Examples include the mobility or 
the thermal conductivity. 

It is well-known that it is possible to approximate 9° (Jjj,j /dvr a) by considering a simulation at 
A = 0. For example, the celebrated Green-Kubo formula [SlIIl] writes (see Theorem 1371 in Section!?] 
for a proof in our specific context): 

f fdnx = r [/(Xo) (VV • dlFx - V • d°xFx) (X°)] ds 
JlRd Jo 

where the subscript ttq indicates that the initial condition Xo is distributed according to ttq. The 
derivative dx (J^d fdT^x) can thus be approximated by considering infinite-time integrals of auto¬ 
correlation functions for the stationary process at A = 0. This formula can be used to approximate 
dx (fgd fdTTx) numerically, which requires at least in some cases to be careful when choosing the 
truncation time in the integral, see for example [7]. Let us also mention another technique discussed 
in |34| . based on the use of Malliavin weights and the Bismut Elworthy Li formula (see 131 dans])- 
In this work, we are interested in so-called non-equilibrium molecular dynamics (NEMD) methods 
which consists in simulating two trajectories with A = 0 and A = £ small, and then considering the 
finite difference when e —>■ 0 (see for example [SlIH])- 


d°x 






when £ —>■ 0 and t —>■ oo. 


Note that the consistency of this estimate is based on the ergodic properties m and To reduce 
the variance of the computation, it is natural to use the same driving Brownian motion for the 
two processes (Xf)s>o and {X°)s>o (see m) and we therefore end up with the natural following 
estimate: ^ 

dx (^J^ fdnx^ ~lj t ->• oo. (6) 

As will be shown below (see Proposition [19]), it is easy to simulate dx{f{Xl)) by using the formula 

d°xif{Xl)) = Tt ■ V/(X°) 

where the so-called tangent vector Tt G R'* is defined by 

Tt = dlxt 


Indeed, the couple {Xt ,Tt) is a Markov process which satisfies the 
stochastic differential equation ©: 

dXf = -XV{X?)dt + V2AWt , 
dTt = {dlFx{X°t) - VW{X°)Tt) 

with initial conditions Xq = Xq and To = 0 (since, we recall, Xq does not depend on A). The 
estimate thus leads to a practical numerical method to evaluate the derivative 9° (/gd fdTJ'x)- 
The main theoretical result of this paper consists in exhibiting sufficient conditions such that 
the following equalities hold true (see Theorem 1401) : 

dx ( f fdTTx) = lim j f d°x{f{X^))ds a.s. (8) 

KJmd ) I Jo 


following extended version of the 


dt, 


(7) 
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and 


aS (^1^^ /dTr^) = ^lim E [dlifiX,^))] . 
Therefore, two natural estimators of 8° (/^d fdivx) are 


-f 

t Jo 


d°x{f{X^))ds 


and 


E a' 




(9) 


( 10 ) 


( 11 ) 


The second estimator is derived from the expected ergodic property on the time marginals: limt-»oo E(/(Xi^)) 
fjgd /dTTA- For both estimators, Theorem |40] can be seen as a rigorous justification of the interversion 
of the derivative dx with the limit limt_>oo and an average in time for the first estimator and over the 
underlying probability space for the second one, since 9° (/gd fdTfx) ~ dx ^limt_>oo j fo f(Xj)ds) 

and dx (fgd fdnx) = 9° (limt-»ooE [/(X^^)]). In addition, we also study the variance of the random 
variable dx(f(Xt)) ~ Tt ■ Xf{X°) which influences the statistical errors associated with the two 
estimators (nni) and GH). 

The proof of ([8]) and ([9]) is based on two main ideas. First, the long-time limit (in law) of the 
couple (Xt,Tt) is identified using a time-reversal argument (see Lemma |42l) in the spirit of the 
argument used in m to study the long-time behavior of two interacting stochastic vortices. We 
are then able to identify the long-time limit of the estimators using the Green-Kubo formula which 
we prove in our setting in Section |4] Second, the justification of the interversion of the derivative 
with the long-time limit and the integrals requires some integrability results, which are based on 
the study of the long-time behaviour of E £qj. _ minSpec(V^F), where (L(f)s>o 

satisfies m with X as an initial condition: 

dYt" = -XV{Yt") dt + x/2dWu 


Xrf = X. 


( 12 ) 


Let us emphasize that we prove all these results in a rather general setting: the state space is non 
compact (namely R'*), the coefficients are only assumed to be locally Lipschitz and the potential V 
is not necessarily strictly convex. The study of the long-time behaviour of the couple {Xt,Tt) is 
very much related to the study of the long-time behaviour of the couple {Y^, DY^), (see Lemma [2lTl 
which may be also useful to analyze other related numerical methods, see [30]. 

The paper is organized as follows. In Section [2| we give preliminary results on the stochastic dif¬ 
ferential equations o and 0, in particular on their ergodic properties and the long-time behaviour 
of the associated Kolmogorov equations. In Section [3| we then introduce the tangent vector Tt and 
study its integrability. In Section [I] we derive and prove finite-time and infinite-time Green-Kubo 
formulae. We are then in position to prove the long-time convergence of the estimators m and m 
in Section [S] Finally, the theoretical results are illustrated through various numerical experiments 
in Section |6| 

In all the following, we assume that Assumptions (Pot) and (Drift) hold, and we do not mention 
them explicitly in the statements of the mathematical results. 


2 Preliminary results on ([I]) and (Ej) and the associated 
Kolmogorov equations 

In this section, we introduce partial differential equations related to the stochastic differential equa¬ 
tions m and m, and study their long-time behaviors. These preliminary results will be crucial to 
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analyze the numerical methods aimed at evaluating (fgd fdrvx) that we study. 

Let us introduce a few notations. For any positive Borel measure /r on R'*, we denote by L^(/r) 
the space of real valued measurable functions on R*^ which are square integrable with respect to p. 
We denote by lLQ(e~^^^^dx) the space of zero-mean square integrable functions: 

Lq = |/ G , J f{x)e~^^^^dx = o| , ( 13 ) 

and by lH[^(e“'^^'“Ma:) the first order Sobolev space associated with the measure e~^^^^dx: 

where V/ is to be understood in the distributional sense. 

2.1 About the solution to ([T]) and the regularity of the law of 

Let us first start by an existence and uniqueness result for the process {Xt)t>o solution to 0. 
Lemma 1. There exists a unique strong solution to ©■ 

Proof. The proof is rather standard. The existence of a weak solution to © is obtained thanks to 
the Girsanov theorem and the existence assumption for the process )t>o (see Assumption (Pot)- 
(iii)). Indeed, 

dXt° = -W(Xt°) dt -f ^/2dWt 

= Fx (X“) dt + V2d(^J^ ■ 

Indeed, under the probability Q such that, for all t > 0, {{Ft)t>o being the natural filtration for 
iWt)t>o), 

SI.,= iri! -\l‘ IF. - F.r(X.«) d,) 

the process Wt = fg -^(-P'o — F\)(X°) ds -I- Wt is a Brownian motion and therefore the triple 

(Xt, VFt, Q) is a weak solution to HI). Note that thanks to the Assumption (Drift)-(i), the Novikov 
conditions are satisfied which justifies the use of the Girsanov theorem [211 Section 3.5-D]. 

Moreover, it is standard to check that trajectorial uniqueness holds for the stochastic differential 
equation o, since from Assumption (Drift)-(ii), x i—>■ F\{x) is locally Lipschitz (see [211 Section 
5.2-B, Theorem 2.5]). As a consequence, by the Yamada-Watanabe theorem (see for example [211 
Section 5.3-D]), the stochastic differential equation ([T|) admits a unique strong solution. □ 

In the sequel, we will need some results about the Radon-Nikodym density of the distribution of 
the process with respect to the equilibrium measure. These properties are given in the two following 
Lemmas. 

Lemma 2. Whatever the choice of Xo, for each A G [0, Ao] and t > 0, Xt admits a positive density 
with respect to the Lebesgue measure on R'^. 

Let us now consider A = 0 and {Yf^)t>o solution to d. For all t > 0, the law of admits a 
density y i—>■ p{t, x, y) with respect to the Lebesgue measure which satisfies the reversibility property: 

yt > 0, e~^^^^p{t, X, y) = U, x),dx ® dy-a.e.. (14) 
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Remark 3. A well-known corollary of d is that if Xq is distributed according to ttq, then the 
process {Xt)t>o solution to Q is reversible: for any t > 0, 

(-^s)sE[o,t] has the same law as 

Proof. Let i/) : —>■ R be a bounded measurable function. By the Girsanov theorem, 

E[i;{X,^)] = E [v^(Xo + FAXo+V2w.)^w.-i |F,(x„+y2w.)|=d.J ^ 

Here, the assumptions of the Girsanov theorem are satisfied. Iirdeed, according to |29i Theorem 2.1] 
these assumptions are satisfied if global-in-time existence and uniqueiress in law hold for both Equa¬ 
tion © (see Lemma [T]) and its driftless counterpart 

dTt = V2dWt, 

which is a mere Brownian motion. The first assertion of Lemma [5] is thus proved. For A = 0, we 
follow |161 page 91] to deduce that 

E]'!/)(X°)] = E ^ii{Xo -b /o*(2Av-|vv|2)(Xo-rv^iv.)dsj ^ 

Now, if one considers the Brownian bridge: 

Vs € (0, t], = x + V2W, + ^(y-x- V2Wt) 

one obtains by conditioning with respect to Wt'. 

E(V'(T")) = E (i}{x -b y2m)e"5'"G+'^^‘)e5'^(“^)c,(a:,a; -b V2Wt)^ 


where g{x,y) =E )(Ss’“)ds^^ This shows that admits a density with respect to 

the Lebesgue measure: 


Pit,x,y) = e 


e ‘u 
(47rt)‘^/2 


From the formula d, it is straightforward to check that 

X, y) = e~^^'"'‘p{t, y, x) 


(16) 


(17) 


by using the fact that g{x,y) = g(y,x) which is a direct consequence of the fact that (Rf’“)sg[o,t] 
has the same law as (Rt(l^)se[o,t]- This concludes the proof of Lemma[2l □ 

Let us now state a few additional results on the dynamics when A = 0 and when (X°)t>o starts 
from a general random variable instead of a deterministic point. 

Lemma 4. Let Xq be distributed according to some probability measure po, o,nd let {X^)t>o evolve 
according to Equation @. Denote by pt the distribution of the random variable X^. 

For all t > 0, pt has a density r{t, ■) with respect to dwo = dx: 

pt{dx) = r{t,x)e~'^^^^ da;. 

Moreover, for 0 < s < t, for da:-a.e. x € R"^, 


r{t,x) =E{r{s,YfL„)) 


(18) 
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where Yf satisfies m- Equation GHl) holds for s = 0 if fj,o has a density r(0, •) with respect to 
duQ = Ax. 

If there exists s > 0 such that |lr(s, •)||l=° < oo, then, for all t > s, 

essinf^gj{dr(s,a:) < essinf^ggdr(t,a;) < esssup^g]gdr(t,a:) < esssup^gRdr(s,i). (19) 

Finally, for any q G [1, oo), if there exists s > 0 such that r{s, •) G L'*(e“'^^“'^da;), then for all t > s, 
r{t,-) G {e~^Ax) and 


yt>S, ■)llL9(e-'^(“)dx) — ■)llL9(e-'^(®)dx)- (20) 

Proof. Let i/) : —>■ R be a bounded measurable function. By conditioning with respect to Xq, and 

usiirg the function p{t,x,y) introduced in Lemma [2] 

^{i’{Xt)) = j J tpiy)p{t,x,y)AyAp,o{x) ( 21 ) 

= J fi’{y)e~''^^'’e^^^'' J p{t,x,y)AiJ.o{x) Ay. 

This shows that the law pt of Xt is r{t,y)e~^^^'’ Ay with 

r{t,y) = J p{t,x,y)Apo{x). 

Likewise, for any s € [0, t], by conditioning with respect to it is easy to check that 

¥,{ip{Xf)) = J 'tp[y)e~'^^^^ f p(t — s,x,y)Apsix) Ay. Now by taking 0 < s < t and using the 

reversibility property (HI, we get 

E(V'(^°)) = y ipiy)e~'^^^'> j e^'^'^^pit-s,x,y)Afis{x)Ay 

= J ilj{y)e~'^^^^ J e^^^^p{t — s,y,x)r{s,x)e~^^^^AxAy. 

Since the law of Xt is r{t,y)e~^^^^ Ay, this shows that, 

Ay-SL.e., r{t,y) = J p{t - s,y,x)r{s,x) Ax = E{r{s,Yt^_J). ( 22 ) 

This integral is well defined since x i—>■ pit — s,y,x) and x i—>■ r{s,x) are non negative measurable 
functions. This shows formula m- 

The maximum principle (I19II is then a direct consequence from this representation formula (I22II . 
Finally, if r(s, •) G E'^{e~^^^^Ax), then, for t > s, r{t, •) G L'^(e“'^^“’^da:) since (using the fact that 
X I—>■ p{t — s,y,x) is a probability density function and the reversibility property (I14II 1 


J \rit,y)\ 


r\,, = 


Ay = J J pit — s,y,x)ris,x)Ax Ay 

< J J pit — s,y,x)\ris,x)\'’Ax Ay 

= J Jpit — s,x,y)e~^^^^\ris,x)\'’AxAy 
= J e~'^^“^^|r(s,a:)|'^d2; < 00. 


□ 


Remark 5. In AvpendixVAl we discuss a stronger assumption on V under which we are able to get 
more precise bounds on pit,x,y). 
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2.2 A Feynman-Kac formula and the Fokker-Planck equation 

For two measurable functions / : M"* —>■ M and (/5 : R'* ^ R, with (/5 locally integrable with respect to 
the Lebesgue measure, consider the Kolmogorov equation associated with the infinitesimal generator 
of the stochastic differential equation ©: 

J dtu{t, x) — Au(t, x) + Fx{x) ■ \/u(t, x) — ip{x)u{t, x), t > 0, x € R"^, 

[ u{0,x) = f{x), a: € R”^. 


In all this section, A is a fixed parameter in the interval [0, Ao]. 

In the following, we will consider solutions to Equation (123) in the following weak sense: 

Definition 6. Let u be a function in the space 

L°° (^[0,T],L^(e"'^<"^da:)) n (^[0, r],L^(|v5(a;)|e"'^^"’d2:)) n (^[0,T],E[^(e“'"^"Ma;)) 

for any T > 0. For f G L?{e~^Ax), we say that u is a weak solution to (1251) if u(0, ■) = f and 
for any function v in (“I (^(p{x)\e~^Ax^, 


_d 

dt 


f u{t,x)v{x)e ^^^^Ax = — f \7u{t,x) ■'Vv{x)e ^^^'^Ax— f (p(x)u(t, x)v(x)e ^'•^^Ax 
+ [ (F\{x) — Fo{x)) ■\7u{t,x)v{x)e~^^^^Ax 


(24) 


in distributional sense. 


Note that the last term in (|24l) is well defined since for A £ [0, Ao], from Assumption (Drift)-(i), 
||Ea — -Fo||Loo(]jii) < CX. Moreover, note that the condition m(0, ■) = f makes sense, since a function u 
satisfying 

u G L2([0,T],tf(e"'^("^da:)) and dtu G L=([0, Tj, IH-i (e-'^^^^da;)) 
actually lies in C{[0,T],L,'^{e~^^^'^Ax)) (see for example [311 Lemma 1.2 p. 261]). 

Proposition 7. Assume f G Ax) and that the function <p is locally integrable with respect 

to the Lebesgue measure and bounded from below. Then, Eguation (123 admits a unique solution in 
the sense of Definition\^ Moreover, this solution is m C([0,+oo),L^(e“'^^“^^da;)). 

In addition, if f G n {^ip{x)\e~'^^^'^A:]^, the solution u is more regular: for 

any T > 0, 


u G L°° ^[0, r],H^(e '^^^^dx) n L^(]v5(a;)je '^*'“^^dx)^ n ^[0, r],L^(e '^^“^^dx)^ . 


Proof. By Assumption (Drift)-(i), there exists Co > 0 such that \\F\ — Co]]L=o(iRa) < Cq. Let C be 

i““i 

a positive constant such that (/? + C-^ is nonnegative. From |31l Lemma 1.2 p. 261], one can 


take e x) as a test function in (I24II and obtain the following estimate 

=—e~^* f I Vw(t, — e~*^* f \u{t,x)\^{ip{x) + C)e~^^^'^dx 

Jlt<l JK'i 

+ e~*^* f {F\ — Fo) ■Vu{t,x)u{t,x)e~^^^'*Ax 

f \Vu{t,x)\‘^e~^^^'^dx — f \u{t,x)\^{ip(x) + C)e~'^^^'^dx 

Jv.'^ J’Sl'^ 


< —e 


2 

„-ct 


< — 


f |Vu(f, a:)|^e ^^^^dx+^—- — Cq f |M(t,a:)|^e 

aR<i 2 Jjjid 

f I Vu(t, a:)|^e~'^^^^dx — f e~‘^^\u{t,x)\^ fip{x) + C — e~'^^^^dx. 

JMd Jurf \ 2 / 


Therefore, by integrating in time, one obtains the following estimate: 



From this estimate, the uniqueness result follows from linearity by taking / = 0 in (I25II . And 
thanks to this a priori estimate, existence can be proved by using a Galerkin method on a countable 
family of smooth functions dense in n L^(|<y9(a;)|e“'^^'^^da:), which exists since the 

measure \ip[x)\e~^^^'^dx is finite on compact sets (see for example |25l chapter II, Theorem 3.5]). 
As explained above, the fact that u G C([0, T], L^(e“'^^“^^da;)) is then a standard result, see for 
example |31l Lemma 1.2 p. 261]. 

In order to obtain the additional regularity, let us take dtu{t, x) as a test function in (I24II : 


[ \dtu{t, 


®)l^e 


-''^^^dx = - 


l_d 
2 dt 


f ]VM(t,a:)]^e f \u{t^x)\^ip{x)e ^^^^dx 

J^d 2 dt jgd 


+ / [Fx — Fq) ■ Vu{t,x)dtu{t,x)e 

JRa 

+ C'o / ]Vu(t, x)]^e“'^^“^^da: + i f ]9tM(t,x)]^e“^^^^dx. 

jRd- 4 J]gd 

Therefore, for a constant Ci > 0 such that ip + Ci is nonnegative, 

/ |Vu(t,x)]^e“'^^“^^dx +/ \u{t,x)f{<p{x) + Ci)e~^^'^^dx+^ f \dtu{t,x)fe~^‘''^^dx 
I at J^d z at J^d 4 J^d 

< Co f ]Vu(t, x)]^ + Cl f u{t,x)dtu{t,x)e~^^^^dx 

dRa JMd- 

< Co f ]Vu(t, x)]^ + Cl f \u(t,x)\^e~^''^^dx + ^ f ]9tu(t, x)]^e~'^^^^dx. 

Jr^ Jr^ 4 Jjjd 
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Using Gronwall’s Lemma, one obtains the estimate after integration in time: 


e f |VM(t, a;)|^e '^^^^da; + e f |u(t, a;)|^((p(a;) + (71)6 

Jm<i ina 

+ f f e~^'^““|9tu(s,a:)|^e~'^^“^^d2;ds < f |V/(a:)|^e~'^^“^^d2; 

Jo JRa JRa 

+ [ |/(a:)|^(¥’(a;) + (7i)e~'^^“^^dx + 2(7i f f e~^‘^°''|u(s,ds. 

Jm‘^ Jo Jm‘^ 

The last term is bounded from above over finite time intervals by a constant times \f{x)\^e~^^^'^dx 
thanks to llSl). Again, this a priori estimate can be made rigorous through a Galerkin procedure, 
and yields the additional regularity stated in the Proposition (see for example [261 Proposition 
11.1.1] for a similar reasoning). □ 

Proposition 8. Let u be a solution to the partial dijjerential equation (1231) in the sense of Defi- 
mtion\^ and assume that the initial condition f is of class with locally Lipschitz second order 
derivatives. If ip is locally Lipschitz, then u, dtu, Vu and \7^u are continuous functions and u is a 
classical solution to (I23j. 


Proof. We use a bootstrap argument based on LjL| regularity results for parabolic partial differen¬ 
tial equations. In order to apply standard results which require 0 as an initial condition, we consider 
V = u — f, which satishes (in the weak sense, see Definition O the partial differential equation: 

J dtvft, x) = Av{t, x) -\- F\(x) ■ Vii(t, x) — p{x)v{t, x) -\- g{x), t > 0, x € R'*, 

\ u(0, x) = 0, X € R"^, 

where 

g{x) = A/(x) -f fA(x) • V/(x) - p{x)f{x) 

is a locally Lipschitz function. 

In this proof, we use the following notation 


n LP([0,T],L«(R‘*)), 

T>0 


where L'^(R‘*) is the L'^ space associated with the Lebesgue measure. We will also use the nota¬ 
tions Lfwhere W stands for the usual Sobolev space. We moreover introduce the notation 


L°°- = Pi ll 

2<(3'<oo 


Last, we set 

= {m e L?L«, dtu £ LfU}. 

Let X be some function in the space of smooth, compactly supported functions on R'^. The 
function xu satisfies, in the weak sense, the equation 


where 


f dt{xv) — A{xv) = il>^ on (0, -boo) x 
1 (xt^)(0,x) = 0, xeR”^, 


< 1 >^ = ixFx - 2Vx) • Vw - (Ax + XP)^ + XQ- 
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From parabolic regularity results, see for example |33l Theorem III.l], one has the implication: 


e L?L^) ^ {xv G n 


( 27 ) 


In addition, from the definition of one has 

(Vx e Co°°, XV e ^ (Vx e c^, e (28) 


since ip and Fx are locally Lipschitz functions (see Assumption (Drift)- (ii))- 

Now, by Definition[6l the function v lies in L°°([0, r],L^(e '^^“^^d®)) n L^([0, T], (e '^^^Ma;)), 
so that <F^ lies in LfL^, for any function x £ Co°- First assume d > 1, and let p be the supremum 
of all those p such that lies in for all x- Assume p < oo. Since 

dp 

—— < mm(p,d), 
d + p 

one can find some p G such that "F^ belongs to for any x (note that > 1, 

since 3 + |<| + | = l)- From (I27II xv lies in L? W^’^, and hence from (I28F iF^ lies in L? 
However, Sobolev embeddings yield 

"F^ G Lt C , where ^ > p since p > ^ 

d — p d + p 


which contradicts the definition of p. As a conclusion, <F^ lies in LjL^“ for any x- In the case d = 1, 
one can directly deduce from G that xv is in Lj C L?for any x, and thus <F^ G 

LtL^“ for any x- In any case, xv lies in L?for any x- 

Now consider the equation satisfied by X^iV, for any coordinate i. One obtains 


dtixdiv) - A{xdiv) = 

(xdiv)(0,x) = 0, 


on (0, -boo) X 
X G M'*, 


(29) 


where 


= (X-P'a - 2Vx) • V(9in) + (x^iFx - (Ax + X‘P)e-i) ■ Vn - {xdiip)v + xdiQ- 

Since xv lies in for any x £ Co° , the function 'F^ is in from the boundedness 

of xdiP, xdig and x^iFx- Then, parabolic regularity (I27II for the heat equation (I29II implies 

Vi, Vx £ Co°“, xdiV £ L?wl’°°- n 

In particular, for any x £ Co°i XV is in From Sobolev embeddings, we deduce that xv 

lies in for any e in (0,1) (C® stand for Holder spaces). From the Holder regularity of the 

initial condition, Holder regularity theory for the heat equation now yields the desired regularity 
on V (and thus on u), see for example |221 Theorem 10.3.3]. □ 

Proposition 9. Assume f G {e~^^^'^dx) and that p is locally Lipschitz and bounded from below. 
Then the solution u{t,x) to the partial differential equation (1231) given by Proposition^ admits the 
following probabilistic representation formula: for all t > 0, 

dx-a.e., u{t,x) = E j^/(y'f'^’“^)e~(30) 
where (Yt^’^)t>o is defined by (notice that {Y(^’^)t>o = {Yt)t>Q is defined by (11211 ) 

r dT/-" = Fx{Y,^’^) dt + V2dWt, 

\ V+ = 
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Proof. Step 1-. Let us first prove a maximum principle for solutions to Equation (I23II in the sense of 
Definition 1^ Assume that the initial condition / of Equation (1231) is bounded from above by some 
nonnegative constant M. Let C be a constant such that tp-\-C is nonnegative. From 1311 Lemma 1.2 
p. 261], and | 171 Lemma 7.6], one can take x) — M)'^ as the test function in the weak 

formulation (1^ . and obtain (using the fact that from Assumption (Drift)-(i), ||Ea—E bUj^oofj-d) < Co 
for some Co > 0) 


f M)+fe 

z at J^d 

= — f \V{e~^*u(t,x) — M)^\'^e~^^^'^dx 

— f {ip(x) + C)e~^^u{t,x)(e~'^^u{t,x) — M)'^e~^^^^dx 

JvA 

+ [ (Fxix) — Fo{x)) ■\'{e~^*u{t,x) — M)'^{e~'^^u{t,x) — M)^e~^^^^dx 
</ |V(e“®'‘w(t,x)-M)+|^e"'^^"’dx 

^ jR'i 

— f {ip{x) + C)e~'^^u{t,x){e~^*^u{t,x) — M)'^e~^^^'’dx 

Jr'^ 

+ x) — M)’^|^e~'^^“’^dx. 

^ jRd- 

By using Gronwall’s Lemma, one therefore obtain after integration in time: 

r 

- \{e-^\{t,x)-M)+fe-^^^'>dx 

JRd- 

^ f f e~®'°'’|V(e~®''’M(s,x) — M)^|^e~'^^^^dxds 

2 Jo Jr'^ 

f f e~'^°‘‘{ip{x) + C)e~'^‘‘u{s,x){e~'^‘‘u{s,x) — M)'^e~^^^^dxds 
Jo Jw^ 

\ [ |(/(x)-M)^ 

^ ./Rii 


2 

< - 


''+l^'=-'^(“^)dx < 0 


so that the function u(t,-) is bounded from above by for any positive t. By a similar 

argument, if / is bounded from below by —M, with M nonnegative, then u{t,-) is bounded from 
below by for any positive t. 

Step 2: Let us now prove the Feynman-Kac formula (|30j under the assumption / G C°°nL°“(R‘^). 
Let X G f > 0 and M > 0. Let tm be the first exit time from B{x, M) (namely the ball centered 
at X and of radius M) for the process (E/’“^)s>o- Since s i—>■ is continuous, tm goes to oo as M 
goes to 00 . Let us consider the solution (f, x) i—>■ u{t,x) to (I23E which is with respect to t and 

with respect to x thanks to Proposition [8l Applying Ito’s formula to u(t — s, V/’^) in the time 
interval [0, t A tm], one obtains 


u{{t - tm)^,yF^ )e 


-/c 


tOTM 


V{YJ 


= u{t, x) + V2 Vu{t - s, Y, 

Jo 




dVE. 


On the interval [0, f A tm], the integrand in the stochastic integral remains bounded, so that this 
integral has zero mean. Taking the expectation, one obtains 


u{t,x) = E |u((f - tmT, 


T(.Yf 


“)ds 


1^ 
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By the above maximum principle the function u is bounded on [0,t] x R'*. With the lower bound 
on ip, the dominated convergence theorem yields, letting M —>■ oo, 

u{t,x) = E _ 

Step 3: Let us now assume that / is in L°°(R‘^). Let f„ be a sequence of C°° functions such 
that sup„>j ||/ii||L°°(R'i) < |l/||L“(Kd), and converging to / almost everywhere. In particular, /„ 
converges to / in by Lebesgue’s theorem. Therefore, the solution Un to Equation (I23II 

starting from fn is such that converges to u(t,-) as n —^ oo in from the a 

priori estimate (1^ . Moreover, one has 


yx € R'*, u„{t,x) = E = E[/,(y,^’")r(n^’")], 

where L is a bounded function satisfying r(yj^’“^) = E ^e~ _ Pqj. remaining of 

the proof, we assume that t > 0 (the formula (I30II clearly holds for t = 0). From Lemma [2] the 
distribution of admits a density p^{t, x, y) with respect to the Lebesgue measure. Therefore, by 
the Lebesgue theorem, E[/„(y/’^)r(y(^’“)] converges to E[/(y/’^)r(yf^’^)] as n ^ oo. This shows 
the equality u{t,x) = E |^/(yj^’“^)e“‘'’Ws ’ )dsj g j._ 

Step 4- Let us now assume that / is in h'^{e~^^^^dx) and let us write f = f'^ — f~ where 
/"*' = max(/, 0) and f~ = max(— /, 0). For n G N, the functions = min(/"'",n) (resp. /T = 
min(/“,n)) are in L°°(R'^) and converge in to /+ (resp. /“). Let us consider the 

solution Un (resp. u~) to Equation (1231) starting from (resp. f~). Since — fn converges 
as n —>■ 00 to / in L^(e~'^^“^^da;), for every t > 0, ut.{t,-) — Un{t,-) converges in L^(e~'^^“^^d 2 ;) to 
u(t,-), where u is the solution to Equation (1231) starting from /. Moreover, one has 


da;-a.e., ut{t,x) = E (y/ 


o-/o ¥’(yr'")d. 


=E[/^(y/’")r(y/’")], 


where r is the bounded function defined above. By the monotone convergence theorem, E[/^(yt^’“)r(yt^’“^)] 
converges to E[/^(y/’“^)r(yj^’“^)]. This shows the equality u{t,x) = E |^/(y/’^)e“for 
dx-a.e. X. □ 


As a corollary of the previous result, we obtain that the law of Xt satisfies a partial differential 
equation (the Fokker-Planck equation). 

Corollary 10. Let Xq be distributed according to some probability measure yo, o,nd let (A°)t>o 
evolve according to Eguation Let us assume that yo has a density ro with respect to dno = 

da: such that ro G L^(e“'^^“^^da:). Denote by yt the distribution of the random variable Xt, 
and by X r{t, x) the density of yt with respect to dno = da; which exists by Lemma^ 

Then, r{t,x) is the unigue solution to the partial differential equation 

J dtr{t,x) = Ar{t,x) — W{x) ■ Vr{t,x), t > 0,x £ R*^, 

\ r{0,x) = ro{x), a: G R'*, 

in the sense of Definition\^ 

Proof. From Lemma [4] we know that 

yt > 0, da:-a.e., r{t, x) = E(ro(yt“^)). 

The conclusion is then a consequence of the Feynman-Kac representation formula (I30j- □ 
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2.3 Long-time behavior of the partial differential equation fl2^ when 

A = 0 

In this section, we are interested in the long-time behavior of the partial differential equation (I23II 
when A = 0, which is related to the stochastic differential equation through the Feynman-Kac 
formula ((301). 


2.3.1 The case ip = 0 

To study the long-time behavior of the solution to (1231) with A = = 0, we introduce the following 

hypothesis (dehned for any 77 > 0). 

Assumption (Poinc(77)). The measure e~'^ satisfies a Poincare inequality with constant rj > 0; 
for any function v in n 


rj f |i!|^(a:)e f |Vu|^(a;)e (33) 

Recall that denotes the functions in with zero mean with respect to 

TTo (see (US])). 

Proposition 11. Let Assumption (Poincf?;)) be satisfied for some positive rj, and let u be a 
solution to (1231) in the sense of Definition O in the case \ = ip = Q, with an initial condition 
f G {e~^Ax). Then u converges exponentially fast to the constant function f{x)e~^^^^Ax 
in the following sense: 


Vt > 0, 

u{t,-) - 1 

^ f{x)e '^^^^da; 

< e”"* 

/- 

f f{x)e '^^^^da; 


J 

Rd 

]L2(e-V(a;)d2;) 


/nd 


Proof. Taking the constant function 1 as the test function in (1241) . one obtains that u{t, x)e~^^^'' Ax 
fgd f{x)e~^^^^Ax for any t > 0. In particular, u{t, •) — f(x)e~^^^^Ax € Lo(e“'^^^^da;). In ad¬ 
dition, from [311 Lemma 1.2 p. 261], one can take u(t,-) — /(a:)e as the test function 

in which yields, using the Poincare inequality. 


f [ /(®)« 

^ m jRii jR'i 


-V{x) 


da; 


5 '^^“^^da: = — 


[ |VM(t,a;)| 

aRrf 


2g-V(x) 


da; 





One concludes from Gronwall’s lemma. 


□ 


This Proposition shows that, under the assumption of Corollarv IlOl lnamelv Apo = r(0, a;)e“'^^“^^da; 
with r(0,-) G L^(e“'^^^Ma;)), the density r{t,x) of with respect to ttq converges exponentially 
fast to 1 if (Poinc(r;)) is satisfied for some positive rj. Actually, the convergence of Apt to 
holds in total variation norm for any initial condition po- 

Corollary 12. Let Assumption (Poincf?;)) be satisfied for some positive rj, and let {Xt)t>o evolve 
according to Equation m- Let us assume that Xq is distributed according to some probability mea¬ 
sure po. Denote by pt the distribution of the random variable Xt, and for all t > 0, denote by q{t,x) 
the density of pt with respect to the Lebesgue measure (which exists according to Lemma\^. Then 

lim ||g(t,-) - e“'^||Li(d^) =0. 

t—^OO ' ' 
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Proof. From Equation dUl), for all t > 0, one has 


dy-a..e.,q{t,y) = J p{t,x,y)y,o{dx). 

Let us fix a positive e. Let us consider to > 0 (to be fixed later on) and q^{to,x) a function in 
which is non-negative, with compact support, such that J q^(to,x) da; = 1 and 

\q{to, x) — q^{to, x)| dx < e. 


/ 


To build such a function q^ito, •), one could for example consider for n large enough 
which indeed converges to q{t,x) in L^(dx) when n —>■ oo. Let us define the function 


q^{t,y) = J p{t-to,x,y)q^{to,x)dx Vt > to, Vy G R"^. 

For t > to, q‘^{t,-) is the density at time t of the process solution to ((2]), with 

distributed according to g®(to,x)dx. Let us now set r^{t,x) = x), the density of 

with respect to e~'^^^^dx. Since r®(to,x) = (to, x) € lJ^(e~^''^'’dx), from Corollary 1101 

(s,x) 1 -^ r®(to -|- s,x) satisfies the following partial differential equation (with unknown r) 

J dtr{s, x) = Ar(s, x) — W(x) • Vr(s, x), s > 0, x G R'*, 

\ r(0,x) = r®(to,x), x G R'^, 


in the sense of Definition (6] In particular, from Proposition 1111 since / r®(fo,x)e ^^^''dx = 1, 

Vt > to, ||r®(f,-) - l|L 2 (e-V(rr)d,„) < - l|lL2(e-'t'(-)dx) 

which is equivalent to 

Vt > to, ||<?®(L-) -e"'^IL2(eV(x)d,^) < e“’’^*"*“>||/(fo,-) - 
By Cauchy-Schwarz inequality, we deduce that Vt > to 

||g®(t, •) - e"'^||Li(dx) = [ |q®(t,x) - X 

JEd I I 

< ||<?®(t,-) - e“''|lL2(,v(x)dx) < e“’’^*"‘“>||/(fo,-) - e“'^||L2(,v(x)dx) 
Moreover, we also have: Vt > to 


lk(i,') -/(i, Olkirdx) = / / (q(to,x)-q'^(to,x))p(t-to,x,y)dx 

aRd aRd 


dy 


< 


/ / 

JRd, J^d 


\q(to,x) - g®(to,x)|p(t - to,x,y) dxdy 


= / |g(fo,x) — g®(fo,x)| dx < e. 
JtL'l' 


We thus obtain: \/t > to, 

||g(t, ■) - e“''||Li(dx) < lk(i,-) - a®(L')llLi(dx) + lk®(L-) - e"'^||Li{dx) 

< e-he“'''*"*“^||g®(to,-) -e"’^IL2(eV(x) dx) 

and the right-hand side is smaller than 2e for t sufficiently large. This concludes the proof. 


□ 
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( 34 ) 


2.3.2 The case ^0 

In this section, we are going to investigate the long-time behavior of the function u defined by 

u{t,x) = E , 

for a generic function ip where, we recall, (ys“)s>o satisfies di). When ip > a for some positive 
constant a, u converges to 0 exponentially fast as t —> oo. We now look for hypotheses on p under 
which this convergence is preserved in the case inf p < 0. 

Notice that by ergodicity, almost surely, limt_>oo \ p(Yf) ds = / ‘pe and therefore, 

the almost sure exponential decay to zero is ensured if p{x)e~'^^^^dx > 0, at any rate in 
p(x)e~^^^^dx). The exponential decay to zero in is more complicated to establish. In 
Proposition m below, we prove this exponential decay under a sufficient condition which contains 
the assumption p(x)e~'^^^^dx > 0. 

When p is bounded from below and locally Lipschitz, from Propositionthe function u defined 
by (1341) is solution (in the sense of Dehnition [6]) to the partial differential equation (1231) with / = 1 
and A = 0. As a consequence, the long-time behavior of (1341) is related to the spectrum of the 
operator A — VP ■ V — p which is self-adjoint in l,^(e~^^^'^dx). 

One way to study this spectrum is to perform the change of variable v{t,x) = 
making Equation (I23J become 

dtv = An -f i (2AP - |VP|^ - Ap) v. 

As a consequence, the long-time behavior of v is characterized by the spectrum of the Schrodigner 
operator A-|-1 (2AP — | VP|^ — Ap '), which can be controlled by the Cwikel-Lieb-Rozenblum bound 
(see for example |10ll24l[28| i. Indeed, for d > 3, this bound states that the number N of nonnegative 
eigenvalues of A + IT satisfies 

N < Ld f max(IT(x),Oj'^^^dx, 

where Ld is some constant independent of W. In particular, if there exists e > 0 such that 

f max Te -I-iAP(x) — i|VP(x)|^ — </p(x),0^ dx < ^, (35) 

jR-i V 4 4 J Ld 

then the spectrum of A + 4 (2AP — | VP|^ — Ap) is included in (—oo, —e). 

There are two main concerns with this approach. First, the constant Ld is unknown, so that the 
criterion is not quantitative. Moreover, by Jensen’s inequality, the exponential convergence to 0 of 

E j^e-^4o 

for (5 > 1 implies the exponential convergence of the function u{t, x) in (1341) . However, in some cases, 
the criterion (1551) may apply to 5p for some J > 1 and not to p. We are going to present another 
criterion which does not present these flaws. 

Proposition 13. Assume that (Poinc(77)) holds for some positive rj, and that 


— oo < inf 1/5 < 0, 



V(x) 


dx > 0, 


and — (inf p) 


Ud‘p‘^{x)e 
(Imd ‘P(x)e-'^(^)dx)^ 


< V- 



V(x) 


dx < OO, 


(36) 

(37) 
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Let E = ip{x)e V = (‘P(®) ~ E)^ e '^^^'^Ax and let u be the weak solution to 

Equation (I23II in the sense of Definitions^ for A = 0, with an initial condition f € {e~^Ax). 

The quantity u^e~^ converges exponentially fast to 0 as t ^ oo: 

3C > 0, Vt > 0, f u^{t,x)e~^^^''Ax <Ce~^*, 

jRd' 

with a rate /3 given by 

( E^ + VA /7 e2 + V\^ , 

+ + + 2— > 0. (38) 

Note that the positivity of the rate (1381) is equivalent to the condition dSTj. Moreover, note that 
the left-hand side in condition (1371) is homogeneous of order 1 in ip, unlike the criterion (1351) obtained 
using the Cwikel-Lieb-Rozenblum bound. As a consequence, if the criterion dSll) applies to S<p for 
some real number <5 > 1, then it applies to p, as expected. 


Proof of Provosition \13[ In this proof, for notational simplicity, we omit the time and space variables 
in the integrals, which are all considered with respect to the Lebesgue measure on S'*. 

From | 31l Lemma 1.2 p. 261], one can take u as the test function in (1241) . obtaining 



On the other hand, taking the constant function 1 as the test function in (1241) . 



Therefore, from the inequality 2ab < + b^, 



17 










By combining (I39II and (SB, one obtains for <5 > 0, 



with 


/cN • r , ,,, d 

ci(d) = 77 + mf (/J - — and C2{S) =-V + 


E 


(43) 


We want to find 5 > 0 such that (I42II ensures exponential convergence to 0 of e~^ as t —>■ + 00 . 

If (5 = 0, one has C 2 ( 0 ) = —77 < 0, so we need 5 > 0. We look for 5 > 0 such that both ci((5) and C2{5) 
are positive. 

From Cl (5) is positive if and only if 


5 77 + inf (p 

^ ^ V 


and C 2 {S) is positive if and only if 

A > ^ 

2E V + E2 ■ 

One conclndes by checking that condition (I37II is necessary and sufficient for the interval ^ v+e^ ’ 
to be nonempty. 

For a given J > 0, Eqnation (I42II gives a convergence rate of 2min(ci((5), From the 

definition of ci and C2 , one can see that ci (<5) is nonincreasing and, under (I37II , is nondecreasing 
in 5. As a consequence, min(ci(5), is maximized for 5ci(5) = C 2 (S). This last equation is 

quadratic, and one can check that its unique positive solution is 


giving the rate IMli- 


□ 


One can see that Equation 1371) is necessary and sufficient for the existence of 5 > 0 such 
that Cl (5) > 0 and C 2 ( 5 ) > 0. One can naturally wonder whether introducing more flexibility in 
the inequalities used in the proof of Proposition HU could lead to a weaker condition. Actually, 
keeping track of the positive term |Vupe“'^ in (13911 . using inequality (13311 in (14011 . and using the 
inequality 2ab < 70 ^ + in HD) leads to the exact same necessary and sufficient condition to 
ensure exponential convergence to 0 . 

Remark 14. One can use the theory of large deviations to prove that the long-time behavior of 
quantities of the form 1341) is necessarily exponential, with a rate given by a variational formula. 

Let {Xt)t>o evolve according to Equation a. According to Donsker-Varadhan’s lemma, the 
random probability measure fit defined by the formula 


1 f* 

h-tiA) = -f j Ix^eAds, 


satisfies a large deviation principle with rate function 

r, ^ // WVjf' 7/3/:R‘*^R, i7 = /(a;)e-'"Wd*, 

I(i^) = < Jmd i I 

[ 00 otherwise, 
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see for example Chapter IV.4]. As a consequenee, in the long-time limit, 

-ilog (e = -^log (e 

converges to the constant a defined by 

a = inf f IV'x/JI {x)e~^^^^dx-\- f (p{x)f{x)e~^^^^dx 
I Jn’l ' ' jR'i 

where the infimum is taken over all probability densities with respect to the measure e~'^''^'’dx, from 
Varadhan’s lemma in large deviations theory. By the change of variables g^ = f, a is equal to 

hid I Vg|^ {x)e~''^^'>dx + cp(x)g^(x)e~''^^^dx 

™ f^d g^(x)e-'^C'>dx 

which is the bottom of the spectrum of the operator —A+VV^-V+ip which is self-adjoint in'L^(e~^^^^dx), 
already discussed at the beginning of this section. 

Let us give two examples where the result from Proposition 1131 applies. 

Example 15. A first example is given by the double-well potential in dimension 1, defined by 

V-/{x) = X* — ')x‘ + Ct, 


where 7 > 0 and = In exp(—a;"^ + 7 ®^) da:). IL^e want to apply Provosition \lA to the case 
where the function ip is a multiple of min SpecV^ Viy(a;) (see Section \8.‘A for a justification of this 
choice for ip). In the present case, this writes p-/^s{x) = 5{12x^ — 27 ), where 5 is the positive 
multiplicative factor. 

Denote by rj-f the optimal Poincare constant associated to the potential V-f. As ^ goes to 0, the 
limit potential x* satisfies a Poincare inequality with constant rjo > 0, owing to its convexity. .4s a 
consequence, the Poincare constants 7]-^ converge to a positive limit. On the other hand, as 7 goes 
to 0 , the quantity 


-(inf 


{ImPiAix)e-'^-’(^'idxY 


goes to 0 , since inf goes to 0 while 


converges to some positive constant. .4s 


a consequence, for any <5 > 0 the inequality dszi is satisfied for 7 smaller than some critical value 
depending on S. Notice that the inequalities ([Mil are satisfied for any 7 > 0 since, for any smooth 
potential E : R 1 -^ R, V"{x)e~^^^^dx > 0 holds from a mere integration by parts. 


Example 16. The second example is given by an identically vanishing potential V{x) = 0, with the 
equation considered on the one-dimensional torus, identified with the segment [0, 27r] with periodic 
boundary conditions. The invariant measure is then the uniform measure on the torus. Consider the 
function p{x) = sin(a:) + a with a > 0. In that cases, the mean value of p is given by a, and p{x) 
is not nonnegative for all values of x as soon as a < 1. 

In that case, the Poincare constant is given by rj = 1 and inf p is given by a—1. As a consequence. 
Equation dSTj writes 


(1-a) 


l + 2a^ 
2 q2 


< 1 


The condition is thus satisfied if a > ao where ao is the unique real root of the equation 


3 1 1 

“ + 2«-2 - 0 , 


given by qo — 0.590. .4s a consequence, for a G (ao, 1), one has exponential convergence of (1 34 II to 
zero while the function p is not uniformly positive. 
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2.4 Existence and uniqueness of an invariant measure tta for ([T]) 

In all this section, we assume that Assumption (Poinc(77)) holds for some positive r;. We would like 
to show that the stochastic differential equation © admits a unique invariant probability measure 
that we denote in the following n\, and to give an explicit formula for this measure. Of course, for 
A = 0, we have 

dvTo = 

and one result of this section is that it is the unique invariant measure for ©• We will use ttq 
as a reference measure to build functional spaces, and to construct the invariant measure tta by 
perturbative arguments, using the crucial assumption on the boundedness of F\ + W = Fx — Fo 
(see Assumption(Drift)-(i)): for A £ [0, Ao], 


\\Fx-Foh^^^,^<CX. 

Let us begin with some notation. We denote hy Cx — Fx ■ V + A the generator of the pro¬ 
cess (X^)t>Q. In particular, £o = —W • V -I- A. Also denote 

Tx^Cx-Co = (Fx + VV)-V (44) 

The space n lHl^(e“'^*'^M 2 :) endowed with the symmetric bilinear form 

{u,v)i-^ f Vu{x) ■ \/v{x)e~^^^'^dx (45) 

Jw’- 

is a Hilbert space by Assumption (Poinc(77)). A consequence of the Riesz theorem is that for 
any u £ L^(e“'^^^Ma;), there exists a unique function v in lJQ{e~^^^^dx) n such that 

Vw G Lo(e~'^^“^^d 2 ;) n H^(e~'^^“^^d 2 ;), f Vv{x) ■ Vw{x)e~^^^'^dx = f u(x)w(x)e~^^^^dx. 

This function is denoted v = —Cq^u since when v is smooth, Vv(x) ■ Vw(x)e~^^^^dx = 

— Jjjd £ov(x)w(x)e~^^^^dx. We denote by V{£o) the domain of £o, defined by 

T>{£o) = |u £ Lo(e~'^^^^da;) n lHl^(e~'^^^^da;), £ov £ Lo(e~'^^^^dx)|. 

For a function u £ L^(e“'^^“^^dx), from the Poincare inequality, one has 

^M|lL2(e-v(x)d,i,) < [ \^i£o^u){x)fe~^^^'>dx=-f {£o\)ix)uix)e~^^'"^dx, 

— II'^O '*^llL2(e-'^(=')da:)ll'*^llL2(e-'^(®)dx) 

which implies 

^ll'^O — ll'*^llL2(e-'^(“^)da;)' (^6) 

In the following, we will use the orthogonal projection operator Ho from dx) onto Lo(e“'^^“^Ma;) 

defined by: 

V/£ L^(e"'^^'^'da;), Ho/= / - / /(a;)e“''^“'^da;. (47) 

Let us now explain formally how we obtain an explicit formula for the invariant measure ttx of ©. 

For any test function ip and since TaI = 0, £xno(<fi) dwA = 0. Thus, by considering / = ToHo(:/?), 
for any test function /, -Ca-Cq ^Uo/dwA = 0 which also writes fjg^d(d -I-TaTq ^IIo)no/dwA = 0 
where I denotes the identity operator. This is equivalent to: for any test function /, J^d{I + 
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Ta-Cq ^no)/ ^^dvTo = Jjgd /dTTo which yields (/+ 7a£o ^^ 0 )* = 1, where * denotes the dual 

operator on the Hilbert space As a consequence, we are naturally led to study the 

operator TxC^^Iio defined from dx) to L,^(e~^^^^dx). The aim of the next Lemma is to 

show rigorously that we can define an invariant measure ttx of o by defining its Radon-Nikodym 
derivative with respect to ttq as (/ + (7A/C[7^no)*)~^l- 

We can now state the result concerning the existence of an invariant measure for ©. 

Lemma 17. Let us assume that Assumption (Poinc(7;)) holds for some positive rj. Then there 
exists Ai G (0, Ao] such that for A G [0, Ai], the dual operator I + (fTxCQ^HQY on the Hilbert 
space lJ^(e~^^^^dx) of the operator I + 7AT(7^no is invertible and has a bounded inverse. 

Let us then introduce, for A G [0, Ai], the function gx G dx) and the associated measure 

tta such that 

dTTA = STAdTTo whcve gx = {I + ^Ho)*)”^! . (48) 

The measure tta is a probability measure which is invariant for the process {Xf')t>o solution to (Hj. 

Proof. Step T. Let us first study the operator 7 aH o. From the boundedness assumption on V’L+ 
Fx = Fx — Fo (see Assumption(Drift)-(i)), the definition of and (I46II . for any u G 

IITaTo = f K^’a + VV){x) ■ X{Co\){x)fe-^^^^dx 

< C\^ [ \V{CQ^u){x)\^e-^^"'>dx 

= —C\^ f {LQ^u){x)u{x)e~^^^'^dx 

2-1 A^ 2 

— C'A ||£o '^\\l^^(s-VMdx)\W\\l2{s-V<.^)dx) — C”— ||w||l2(j,-V(x)j2,). 

As a consequence, the operator TxCq^ is bounded from lJ^{e~'^^^'^dx) to with: 

1 / (J 

ll£(Lg(e-''(“=)d3;).L2(e-'A(i«)da;)) — Y (^9) 

By composition, 7 a£q ^Ho is thus a bounded operator from L^(e~'^^“^^da:) to itself, with a norm 
of order 0{\), and so is (TaTi^^Ho)*. As a consequence, for A small enough, the operator I + 
(TaTiC^Ho)* is invertible from L,^(e~^^^^dx) to itself. 

Step 2-. Let us now introduce the function gx G dx) defined by 

<7A = (/ + (rA£o^no)*)-4 

and let us prove that dvTA = flAdvro is invariant for the stochastic differential equation ©■ Let 
be the solution to m with initial condition = x (see m)- Using the Markov 

property, the aim is to prove that, for any C°° bounded test function / : R'^ —>■ R, 

[ ^{fiY/""'))9\{x)e~^^'^''dx= f f{x)gxix)e~^^'^'’dx. (50) 

JlRd JRd 

From Proposition m we know that u{t,x) = E(/(y/’“^)) is the solution to (l23ll (with yj = 0), and 
from Proposition [3 we have for any T > 0, 

M G L°° (^[0,T],e^(e“'''">da:)) n (^[0,T],L^(e"'^^"^da;)) . 
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Moreover, from Proposition [8l m is a classical solution to (1231) . Therefore, 


d_ 

dt 




[ E(/(y/’'")) 5 A(x)e ^ / u{t,x)gx{. 

jR'i m J]jd 

= f dtu{t,x)gx{x)e~^^^^dx 

= [ jCxu{t,x)gx{x)e~^^'^^dx, 

Jv.d- 


and £xu = dtU € L^(e ^^^^dx). Now, notice that for any function i/) which is the sum of a C°° 
function with compact support and a constant, 

Tx^ = TxCq^^oCo^ 


holds true since 7a sends constant functions to 0. Therefore, for any such function -0, one has 




CxtpgxdiTo = f [{Co+ Tx)ip]gxd-Ko = f [{I + TxCo^'no)Co^] gxd-Ko 

Jr'^ JRii 


= [ [^oV’] 

Jm.'i 

Jr<‘ 


{I + {TxCo^nor)gxdno 


Co'tpd'KQ. 


Since ttq is invariant for the dynamics © with infinitesimal generator Co, the right-hand side is zero. 
By density, the equality Cxipgxdno — 0 holds for any function tp such that Cxip G l-?{e~'^''^^dx). 
Therefore, ^ E{f{Y^^’^))gx{x)e~^^^^dx = 0 which yields (l50ll after integration in time over [0, t]. 

Step S'. Let us finally check that tta is a probability measure. First, one has 


f gxdTvo^ f {I + {TxCo^TloTy^l{I + {TxCo^no))ldT:o = f dTro 

jR'i Jk<‘ 


= 1 . 


Second, one can prove that (?a > 0. Indeed, from (fSCTI) and the fact that admits a density 
p^{t,x,y) with respect to the Lebesgue measure (see Lemma [2]|, we have 

[ [ f{y)p^{t,x,y)dygx{x)e~'^^^'>dx = f f{x)gx{x)e~'^‘-'^'‘dx. 

jR‘i JR'i JR‘i 

This equality holds for any smooth test function / and, by a density argument, one can apply it 
to the bounded function f{x) = sgn{gx{x)), where sgn(y) = lj,>o — lj,<o denotes the sign function. 
One thus obtains 

[ [ {sgn{gx{y))sgn{gx{x)) - l)p^{t,x,y)dy\gx\{x)e~^^'^'>dx = 0. 

Jr^ Jr^ 

Thus, (sgn(srA(y))sgn(5iA(a:)) - l)p^{t,x,y)\gx\{x) = 0 da: (g) dy-a.e.. Since p^{t,x,y) > 0 da; (g) dy- 
a.e. (see Lemma O and |pa|(®) da: > 0, this implies that dy-a.e., sgn{gx{y)) = 1 or dy-a.e 
sgn{gx{y)) = —1. The conclusion then follows from the fact that gAdTro = 1. □ 

Notice that in the case A = 0, we indeed have = 1- The next result states the uniqueness of 
the invariant measure for ©• 
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Lemma 18. Let us assume that Assumption (Poinc(77)) holds for some positive rj. For X G [0, Ai] 
("Ai being the constant introduced in Lemma m the unique invariant measure of the stochastic 
differential equation © is the probability measure tta defined by (HSli . This probability measure is 
equivalent to the Lebesgue measure on R'^ and for any initial condition Xo, 

V/ G Li(7rx), P f lim i f f{X^)ds = f fdiTx] = 1. (51) 

Proof. Let A G [0, Ai]. Lemma [T7| ensures that n\ defined by (I48II is an invariant probability measure 
for dXt^ = F\{X}')dt + y/2dWt. For Xo distributed according to any invariant probability measure, 
Lemma [2] ensures that this measure is equivalent to the Lebesgue measure. As a consequence, all 
the invariant probability measures are equivalent and the dynamics admits exactly one invariant 
probability measure nx- Since tta is the only invariant probability measure, it is ergodic (see for 
example |271 Theorem 3.8 and Equation (52)]) and denoting by {Y^'^)t>o the solution to © started 
from Yo = X £ R'*, 

V/ G L^(7rx), da: a.e., V ( Yua \ f /(T/’“')ds = f /dvrx^ = 1. 

For any initial condition Xo, since the law of Xf is absolutely continuous with respect to the 
Lebesgue measure (see Lemma [2)), (I51II follows by the Markov property. □ 


3 Tangent vector of the diffusion 

In all this Section, {Xt)t>o denotes the process solution to ©, with an initial condition Xo which, 
we recall, does not depend on A. We establish various results on the tangent vector Tt defined 
by dSll), which naturally appears in the estimators and m to evaluate 8° (J^d fdnx). 

3.1 Definition and interpretations of the tangent vector 

If the function / is differentiable, one can write 9° = TfV f{Xt), where the process {Tt)t>o 

is the so-called tangent vector, defined as 

Tt = dlX^, (52) 

and the existence of which is ensured by the following proposition. 

Proposition 19. For any t > 0, the function A i—^ Xt is almost surely differentiable, and the 
definition of the tangent vector (EH) makes sense. Moreover, {Tt)t>o almost surely satisfies the 
following ordinary differential equation whose coefficients depend on (A’)’)t>o.’ 

f^=aSTx(A°)-VV(A°)T„ 

I ro = o. 

Proof. By (Drift)-(i) and the continuity of VF and (X°)t>o, t |9x^A(A't')| + |V^F(X°)| is 
locally bounded. Hence (1531) admits a unique solution {Tt)t>o by the Cauchy-Lipschitz theorem. 
Let us prove that for i > 0, A i—>■ X^ is differentiable at A = 0 with derivative equal to Tj . For 
A G [0,Ao], we set rx = inf{t > 0 : \Xi \ > sup^gjQ |X°| -|- 1} with convention inf0 = -boo. Let 
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Lj° = |jfO|+i |VV(a;)|. For t € [0,t\, one has 

ptAr^ . . I I 

sup < / \Fx{X^) + VF(X,^) + VF(X°) - vy(x,^) ds 

se(o,t] io ' ' ' ' 

< C\t + if f Ids 

Jo 


vO 

Cle^t ‘-1 


SO that sup^g[Q j] |X^at>, — -^sataI < —3fD For - 1 one deduces that tx > t 


yO 

C[e^i ^- 1 ) 


and supgg[Q |X^ — | < ^ - -X. In particular, Xt converges to Xt uniformly for t £ [0, f|. 

Now, for t > 0, 

Xt^-X°= fiFxiX^) - Fo{X^))ds + r V^V{i^){X° - X^)ds, 

Jo Jo 

where, by a slight abuse of notations, stands for the matrix {dij and Vi € 


{1,. .., d}, £ [X°, Xi-]. For s G [0, t\ and A G ( 0, Ao A —^ 

Hence for t £ [0, f]. 


C{e^t *-l)J 


, |^s’‘| < SUPig[o_J] \^t I + 1- 


sup 

se[o,t] 


-y-A -yO 
As — As _ 0 

J- .<? 


A 


< 


f 

Jo 


Fx{X^) MXll _ gOp^^j^O^ 


+ Lt 


A 

yA yO 

~ rr^O 

> .s 


ds. 


By (Drift)-(i)-(ii) and the uniform convergence of X^ to X^ for t £ [0, t], 


hm r ^-(^-)-^°(^-) -agFA(X°) 
A^oio ^ 


+ |VV(X°) - VV(C.^)||T°|ds = 0. 


With Gronwall’s lemma, one concludes that sup^gjQ ^ — F? | converges to 0 as A —>■ 0. □ 

We have the following expression of {Tt)t>o as an integral: 

Proposition 20. Define the resolvent {Rxo{s,t))s,t>o associated with Equation ([53ll as the solution, 
with values in to the following ordinary differential equation: 


dtRxo{s,t) = -XW{X?)Rxois,t), s,t > 0, 
Fx° (®i ^) — F, s > 0, 


(54) 


where Id is the d x d identity matrix. The resolvent satisfies the following semigroup property 

Vr, s,i G [0, oo), Rxo{s,t)Rxo{r,s) = Rxo{r,t). (55) 

One can recover the tangent vector from the resolvent through the following formula: 


Vt>0, Tt= r Rxo{s,t)d°xFx{X°)ds. 
Jo 


(56) 
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Proof. The semigroup property (1551) is a consequence of uniqueness for Equation (1541) . satisfied by 
the two processes {Rxo{s,t))t>o and {Rxo{r,t)Rxo{r, s)~^)t>o- 

In view of the differential equations satisfied by {Tt)t>o and {Rxo{s,t))t>o, one has, from the 
equality Rxo(t,0) = Rxo{0,t)~^, 

dt {Rxo (t, 0)Tt) = —Rxo {t, 0)dt {Rxo (0, t))Rxo (t, 0)Tt + Rxo (^i 0)dtTt 
= Rxo it, 0)V V(X°)7?;,o(0, t)Rxo it, Q)Tt 
+ Rxoit,Q)dlF^iX°) - Rxoit,0)X^ViX°)Tt 
= Rxoit,0)dlF^iX?). 

Integrating over [0,t], one obtains 


Rxoit,0)Tt = r Rxo is, 0)dlFxiX°)ds, 

Jo 


and the result follows by using the semigroup property (ESI). 


□ 


Notice that the resolvent is also the differential of the trajectory with respect to its initial 
condition. 

Lemma 21. Let solve m- Then for any t > 0, x is on R'* with Jacobian matrix 

iDYr)ij = given by DY^ = iiv»(0,t). 

Proof. By standard results on ordinary differential equations, x i— Yf^ is with Jacobian matrix 
DYt solving the equation 

rt 


Vt > 0, DYt" = Id- XYiYf)DYfds, 


/ 


obtained by spatial derivation of Y" = x — V’E(yj“)ds + ^/2Wt. By uniqueness for (1541) . one has 
DYt" = RY^iQ,t). □ 

In the following, we will need the following result about the link between the forward resolvent 
and its backward counterpart. 

Lemma 22. Let (41 ^satisfy Eguation dSJ with Yo distributed according to ttq. From the 
reversibility of the dynamics m, the process iZs)o<s<t defined by Zs = Yt-s has the same law as 
iYs)o<s<t, and one has the relation 

RYiO,s) = Rzit - s,t), 

where Ltz is the transposed matrix of the resolvent associated with Z. 

Proof. Uniqueness holds for the ordinary differential equation satisfied by s i—>■ (0, s): 

r dR, 


^^-is) ^ -V^ViY)Ris), 
RiO) = 7d. 


(57) 


One can check that s i— R%it — s,t) also solves (1571) . Indeed, since, by the semigroup property, 
Rzit — s,t) = Rzit,t — s)~^, one has, for s € [0,t], 

dsRzit — s,t) = —Rzit — s,t) idaRzit,t — s)) Rzit — s,t) 

= —Rzit — s,t)V'^V iZt-s)Rzit,t — s)Rzit — s,t) 

= -Rzit-s,t)VWiZt-s) 

= -Rzit-s,t)vWiYs). 
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This concludes the proof. 


□ 


3.2 Almost sure boundedness of Rxo{0,t) and Tt 

The tangent vector can take large values, since the second term in the right-hand side of (|5^ will 
provide exponential growth for {Tt)t>o, typically when the trajectory is close to a local 

maximum of V, or when it crosses a saddle point of V. In the sequel, we need some assumptions 
on V to control this behavior. 

3.2.1 Local-in-time boundedness of Rxo{s,t) and Tt 

Let us first introduce an assumption which will be sufficient to get the local-in-time boundedness 
of Rxo{s,t) and Tt. 

Assumption (min Spec). The matrix-valued function V'^V : R'* ^ is bounded from below, 

in the sense that there exists q G R such that, for all x,h € R"^, 

h-\/^V{x)h > a\hf. 

Equivalently, the spectrum of V^V{x) is bounded from below by a, uniformly in x. 

Under Assumption (min Spec), the random variables Tt and Rxo{s,t) are bounded: 

Lemma 23. One has 

VO < s < t, |Rxo(s,t)| < e-J'>-Spec(vV(x2))d„^ 

being endowed with the matricial norm associated with the Euclidean norm on In addition, 
if the Assumption (min Spec) is satisfied, for any T > 0, the random variables suPq<j,<j<j, \Rxo (S) t)l 
and supQ<(<y \Tt\ lie in L°“(fl). 

Proof. For any vector x, one has 

dt\Rxo{s,t)xf = -2iRxo{s,t)xfV^V{X?)iRxois,t)x) 

< —2 min Spec (V^U(X°)) \Rxo{s,t)x\^. 

As a consequence, one has the estimation 

\Rxo{s,t)xf < |j;|2e-2/^™Spec(v=V(x2))d„ 

so that (I58II holds. If the Assumption (min Spec) is satisfied, this inequality proves that R{s, t) is 
in L°°(fl) locally uniformly in time. From the expression (1561) of Tt and the boundedness of d^Fx, 
one also concludes that Tt G L°°(fl), locally uniformly in time. □ 

3.2.2 Global-in-time boundedness of Rxo{0,t) 

We need some additional assumption on the convexity of the potential for (Rxo(0,t))t>o to be 
bounded globally in time. 

Assumption (Conv). The potential V is such that 

f max (O, — min Spec (V^U(a:))) < oo and f min Spec (V^U(a:)) > 0. (59) 
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In this assumption, the first inequality, always satisfied under (min Spec), ensures that the 
integral min Spec is well defined in (— 00 , 00 ]. We refer to Appendix m for a 

discussion of this Assumption. 

Lemma 24. Under Assumptions (min Spec) and (Conv), the resolvent matrix Rxo{0,t) almost 
surely converges to 0 as t goes to infinity, with exponential rate. Namely, for any ft with 

0 </3 < f minSpec (V^y(a;)) 

there exists an almost surely finite random variable C > 0 such that 

Vt > 0, |i?jfo(0,t)| < Ce~^\ (60) 

Proof. From Lemma [TSl by ergodicity (see dSlIl '). one has 

lim - [ min Spec (V^F(As)) ds = [ min Spec (V^F(a;)) a.s. 

t^oo t Jo JMd 

We conclude by combining this limit with Equation (1581) . □ 


Remark 25 (On the Assumption (Conv)). While Assumption (Conv) is automatically satisfied 
in dimension 1 from a mere integration by parts, this is not the ease in higher dimension. Indeed, 
if one applies the integration by parts formula in this case, one only obtains that 




VV{x)(SVV{x)e~'^’-"^dx 


is a positive definite matrix (because of the integrability of e~^, for any y in the function x 1 -^ 
W{x) ■ y cannot be the zero function), so that the minimum of its spectrum is positive. A coun¬ 
terexample to Assumption (Conv) is given by a tensor potential V{x) = U(xi) + ... + U{xd) with 
a well chosen function U. Indeed, in this case the left hand side of equation dSH) rewrites 


f min (U''(xi))e ... da;^ = E 

min U" (Xi) 




where Xi are i.i.d random variables with distribution e ’^^^'^dx. IfU is chosen so that U" is bounded 


and has a strictly negative minimum, then the sequence 


min U” iXi) 

ie{i.....d} 


converges almost 


surely as d goes to infinity to the negative constant min U". Then, from the dominated convergence 

theorem, the quantity E min U” (Xi) converges to min U", and is thus negative when d is large 
[ie{l,...,d} J 

enough. 

Remark 26 (On the assumptions of Lemma [24ll . Assumption (Conv) is not necessary for (1601) to 
hold. Indeed, if the matrices X^V{x) commute, the matrix Rxo is given by 


RxoiO,t) = e-^o^"nx))as 

and the convergence of i V^V(As)ds to the positive definite matrix V^( 2 :)e“'^^^^da: implies 

that (1601) holds for (3 < min Spec ^ X^V(x)e~^^^'^dx^, even in the cases when V does not satisfy 

Assumption (Conv). An example where the matrices W^{x) commute is the case of a tensor 
potential V{x) = U(xi) + . .. + U(xd). As seen before, U and d can be chosen such that V does not 
satisfy Assumption (Conv). 
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However, it is likely that Lemma \E^ does not hold under the sole ergodicity property: 
lim - f (vV(X°))ds= / VV{x)iS,VV{x)e~^^"'>dx a.s. 

t-»oo t Jo 

Indeed, there exists some family of symmetric matrices (j4t)t>o converging in the Cesdro sense to a 
positive-definite matrix, for which the solution to ^Rt = —AtRt, Ro = Id does not converge to 0 
as t goes to infinity. An example of this phenomenon is given by 


At — 0,t 




where Qt = 


f cost 
\^sin t 


— sint 
cost 


Indeed, the family {At)t>o converges in the Cesdro sense to Id as t goes to infinity, but the associated 
matrix {Rt)t>o diverges. To show this last point, consider the matrix Mt = Rt. Since dtUt ~ 



holds, 
are — 


dtMt = 


1 1 
-1 -3 


Mt 


As a consequence, Rt = fit exp ( t 


1 1 
-1 -3 


. The eigenvalues of the matrix 


1 

-1 


1 — %/3 and — 1 + %/3, the latter being positive, so that Rt diverges as t goes to infinity. 



3.3 Boundedness of moments of Rxo{s,t) and Tt 

In the sequel, we will need to control the moments of Tt. From Equation (1561) and the boundedness 
of dxF\ (see Assumption (Drift)-(i)), this boils down to estimating the moments of Rxo{s,t). For 
this purpose, from (l58l) . it is enough control expectations of the form E , 

where /I is a positive constant. 


3.3.1 Preliminary result when Xq ~ ttq 

One can deduce from Proposition [13] a criterion for exponential convergence of the moments of 
Rxo{ 0 ,t) to 0 as t —>■ oo. To state the result, we need to strengthen the assumptions (min Spec) 
and (Conv) which is the point of the following assumption. For any p > 0, let us consider: 

Assumption (Spec(p)). Assume that 


— 00 < inf minSpec(V^F(a;)) < 0, 


f minSpec(V^F(a:))e ^^^^dx>0 and f (min Spec(V^F(a;)))^ e < oo, 

Jmd Jud 


— (inf minSpec(V^P(a:))) 


Jjjrf (minSpec(V^F(a;)))^e ^^^^dx 
(Jjjrf minSpec(V^F(a:))e“'^(‘^lda:)^ 


< P- 


Notice that for p > 0 and fi > 0 under assumptions (Poinc(p)) and (Spec(p//3)) then the 
assumptions (13611 and (|37|) of Proposition m are satisfied with (p{x) = /9 min Spec(V^P(a:)). 

We are now in position to state a simple consequence of Proposition 

Proposition 27. Let {X?)t>o solve © starting from Xq distributed according to ttq. Assume 
that (Poinc(p)) and (Spec(p//3)) hold for some rj > 0 and /3 > 0. Then there is a constant 
C € (0, +oo) such that 

yt > 0, E,J|i?^o(0,t)|'’] < 
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Here and in the following, the notation means that the initial condition Xo of the processes 
{Xt )t>o solution to © is distributed according to -kq. 

Proof of Proposition \27\ To apply Proposition to the function (p{x) — /? min Spec(V^P(a;)), we 
need this function to be locally Lipschitz. Since is locally Lipschitz (see Assumption (Pot)- 

(i)), this is a consequence of the Lemma [28l given below. By Proposition the function u{t,x) = 
E ^e~ ^mmSpec(v v{Y^ ))dsj jg ggiytion to Equation (1131) in the sense of Definition [S] for ‘fiix) = 
/JminSpec(V^E(a;)) and f{x) = 1. 

Since conditions ([35)1 and (133 hold for this choice of ip, Equation (I58II and Proposition 1131 give 

E-n-o fl^xo(0,i)|^l < [ u{t,x)e~^^^^dx < f f M^(t,a:)e“'^^^^da:^ < 

L J Jmd VJgd J 

for some positive constant C. □ 

Lemma 28. The function A i— ^ min Spec(A) is a Lipschitz function on the space of symmetric dxd 
matrices. 

Proof. Let A be a symmetric matrix, and let a: be a vector in R'* such that |a;| = 1 and min Spec(A) = 
x ■ Ax. Then, for any symmetric matrix B, one has 

min Spec(i3) < x ■ Bx = x ■ {B — A)x + x ■ Ax < | A — i3| + min Spec(A), 

jjdxd endowed with the matricial norm associated with the Euclidean norm on By 

exchanging A and B in the previous inequality, one obtains 

I minSpec(A) — minSpec(i3)| < |A — B\. 


□ 


3.3.2 Uniform-in-time boundedness of moments of Tt 

Numerically, the computation of © through the long-time limit of a Monte Carlo approximation 
of the expression E[rt • V/(X)’)] is only possible if Tt has a bounded variance uniformly in time. 

A case where this fact is easily proved is when the function V is T;-convex, where is a positive 
constant. We recall that this means that the spectrum of X^V{x) is bounded from below by rj, 
independently of x. More precisely, one has the following proposition. 

Proposition 29. Assume that the V is rj-convex, for a positive constant rj. Then, for any a > 1, 

supE|rt|“ < 00. 

t>0 


In particular, Tt has a bounded variance uniformly in time. 
Proof. By (I58II and the boundedness of dxFx, one has 


E[|Tt|“] =E 


r Rxo{s,t)dlFx{X°)ds 
Jo 


a- 

<E 

(c re-“'*-“^ds) 



\ Jo J . 


< 00 , 


so that Tt has a finite moment of order a. 


□ 


The convexity assumption on the potential can be loosened, as shown in the next Proposition. 
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Proposition 30. Let a G [l,+oo). Assume that (Poinc(77)) holds for some positive rj, that the 
initial condition Xq is distributed according to a measure po having a density with respect to the 
measure which is in for some p G (l,oo], and Assumption (Spec(77(p — 

1)/(q:p))) holds (with the convention 77(00 — l)/(aoo) = rilo.). Then, 

supE|rt|“ < 00. 

t>0 


and, when a>2,Tt has a bounded variance uniformly in time. 


Proof. By (I56II and (Drift)-(7), 

E'/“[|rtn < rE^/'^[\Rxo{s,t)dlFx{X°)nds<C f E^/'^[\Rxois,t)nds 

Jo Jo 

Let Pa denote the law of X® for s > 0 and {Yt)t>o be a solution to Q with Yq distributed 
according to tto. We notice that the Markov property gives: for 0 < s < t, 


E[|i?^o(s,t)n =E 


\RY{s,tr 


dps 

Q-Vix)dx 



Using Holder inequality with q — p/{p — 1) {q — 1 if p = 00 ), Lemma [3] and Proposition 1271 one 
deduces that for t > s > 0, 


E[ii?xo(s,t)n<E[|i?v(s,t)in’ 


dps 

e-^(“:)da; 


LP(e 




This concludes the proof. 


□ 


We are now in position to give sufficient conditions for the finiteness of the variance of the two 
estimators (HSl) and mi- 

corollary 31. Let / : R'* —> R 6e a function such that V/ is bounded. Let us assume that 
either V is g-convex (for a positive constant rj), or that there exists rj > 0 and p G (l,oo] such 
that (Poinc(p)) holds, Xq is distributed according to a measure po having a density in L^(e“'^^“^^da;) 
with respect to and Assumption (Spec(p(p — l)/(2p))) holds. Then, 

supVar(Tt • Vf{X°)) < 00 and supVar (— [ Ts ■ V/(Xa)ds^ < 00 . 

t>o t>o V i lo / 

Proof. These results are simple consequences of the boundedness of Vf and Proposition 1301 for 
0 = 2 . □ 

From the Central Limit Theorem for trajectorial averages (see for example |12l Section 2.1.3, 
Theorem 6.3.20]), it is expected that the variance of j JqTs ■ Vf(Xs)ds actually scales like 1/t 
in the limit t —>■ 00 . This requires for example to prove the existence of a solution to the Poisson 
problem associated with the Markov process (Xs,rs)s>o, which does not seem to be ensured under 
our set of assumptions. We leave the study of this issue to a future work. 

Remark 32. Under the additional assumption (V) given in Appendix\^ it is possible to extend the 
previous results to more general initial conditions. Assume that the initial condition Xq is distributed 
according to a measure po such that the measure e^'^dpo can be written as 

e^^^^^djio = fix)dx + dn, (61) 
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where f is some function in L^(da;) withp G [1,2] and v is some finite measure on From (1821) . for 
any t > 0, ytt is absolutely continuous with respect to with dx). Now, 

by the semi-group property satisfied by Rxo, (1581) and the fact that — min Spec(V^y(a:)) < C < oo, 
one has for e > 0, 

|7?xo(s, 1)1 < |l?jco(s V £, t)Rxo{s, s V e)] < |-Rxo('S V e, 1)|. 

For a > 0, using a similar change of measure as in the previous proof, the fact that £ 

and Proposition^^ one deduces that under Assumptions (Poinc(77)) and {Spec{p/{2a))), 

for t > sV £, 


E[|i?jfo(s,l)r] < e^(®"'’>’^E[|i?y(sVe,t)|^“]^E 


dpe 


e-^(Ndx 


{Ye) 


This estimation remains valid for 0<s<t<e up to increasing C, since then, by (|58l) and 
the fact that — mmSpec(V^P(a:)) < C < oo, (s, t)| < In conclusion, for a > 1, under 
Assumptions (V), (Poinc(r;)) and {Spec{ri/{2a))), sup^^Q E[|rt|“] < oo if po satisfies (1611) . 


4 The Green-Kubo formulae 

A first way to compute the derivative m is to use the Green-Kubo formula (see for example |18| 
for a mathematical approach and [6] [14] for physical motivations). This formula gives an expression 
of (O in terms of the time autocorrelations of (A°)t>o, where (X)')t>o satishes (|5]) with an initial 
condition Ao being distributed according to the equilibrium measure ttq. 

4.1 Finite time Green-Kubo formula 

We start with the Green-Kubo formula in Hnite time, which will not be used in the sequel of the 
paper, but motivates the infinite horizon Green-Kubo formula. 

Theorem 33. Let f G lJ^{e~'^^^^dx) be a Lipschitz function and let V/ be its gradient in the 
sense of distributions which can be identified with its almost everywhere gradient. Suppose that the 
initial condition Xq is distributed according to the equilibrium measure no and that Assumption (min 
Spec) is satisfied. Then, for any 1 > 0, for any A G [0, Ao], f{Xt) is integrable and A i— >■ E.,ro [f{Xt)] 
is differentiable at 0 with derivative 

dlE^,[f{X^)] = [v/(Ao) • Rlo{0,s)dlF>,{X°)'\ ds. (62) 

Proof. Since A° is distributed according to e~^^^'^dx, Proposition [12] and the chain rule ensure that 
A I—>■ f{X{') is a.s. differentiable at A = 0 with derivative Xf{X^).Tt. 

To justify the interchange between the derivation and the expectation, we need some integrability 
property. For A G (0, Aq] and 1 > 0, one has, using (min Spec) and (Drift)-(i) for the inequality: 

\X,^ - A°|2 = 2 [\fx{X^) + XV{X^)).{X^ - X°)ds + 2 f {XV{X°) - XV{X^)).{X^ - A°)ds 
Jo Jo 

< C\^t -f (1 - 2a) [ jxi' - A°|^ds. 

Jo 
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As a consequence, 


VAg (0,Ao], 


\X^-X^\ 


< c 


,(l-2a)t _ ^ 


(63) 


A2 - “ 1 - 2a 

with the convention that the last ratio is equal to t if 1 — 2a = 0. With the Lipschitz continuity 


of /, one deduces that the random variable 




is bounded by a deterministic constant not 


depending on A. The integrability of f(Xt) then follows from the integrability of f{Xt) where Xt 
is distributed according to and / G 

Moreover, by Lebesgue’s theorem, A i—E^g[/(X/')] is differentiable at 0 with derivative 


dlE^olfiXh] = E^„[V/(A°) • Tt] = [V/(X°) • ds, (64) 

Jo 

where we used (1561) for the second equality. All terms in Equation (1641) are well defined thanks to 
Lemma [23l Let us now rewrite the right-hand side of (1641) . By introducing the process (Ts)o<s<t = 
(A')’_s)o<s<t (which has the same law as (X°)o<a<t), using a change of variable s —>■ t — s and 
Lemma m we get 

/‘e,„ [VfiX?) ■ Rxois,t)d'iFx{X°)] ds = /‘e,„ [XfiX?) ■ R^o{t - sR)dlFx{Xts)\ ds 
Jo Jo 

= ^‘e,„ [v/(yo) • R'^iO,s)d°xFxiYs)] ds 
= [v/(Ao) • i?Jo( 0 ,s)aSf’A(X°)] ds. 

This completes the proof of dMl)- □ 

Remark 34. The conclusion of Theorem \33\ still holds if f £ l-}{e~^^^'^dx) is a function such 
that V/ is uniformly continuous on and V/ G h^(e~^^^'^dx). 

It is possible to give another expression of the right-hand side in (1621) . 

Proposition 35. Let f G h^{e~^^^''dx) be a Lipschitz function. Assume (min Spec)and consider 
the process (X°)t>o satisfying © with an initial condition Xo being distributed according to the 
eguilibrium measure ttq. For almost every s > 0 

E,„ [v/(Ao) • Rlo{Q, s)aSJ’A(A°)] = E,„ [/(Ao) (VE • d^Fx - V • d°xFx) (A°)] . (65) 


Proof. Since / G L^(e '^''^^dx), by Proposition [7] the partial differential equation 


j dtu{t,x) = Au{t,x) — XV{x) ■ Xu{t,x), t > 0, a; G R, 
u(0,a:) = f{x), a: G R. 

admits a unique solution u in the sense of Definition Moreover, according to Proposition [51 


Vs > 0, da; a.e. , u(s,x) = E[/(yif)], 

where (Yi )t>o solves m- When s > 0, from Lemmas 151 121l andl23l and Assumption (min Spec), 
one can apply the dominated convergence theorem to differentiate E[/(Kf)] with respect to x, 
obtaining V,,E[/(Kf)] = E [R^.{0, s)X f{Yf)]. Since u G nT>o (^[0, r],E[i(e-'^f"Ma;)) , ds a.e., 
m(s, .) admits a distributional gradient denoted by Vu(s,.) and 

ds a.e., da: a.e., Vu{s,x) =E (0, s)V/('Kf )j . 
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When Xo is distributed according to tto, from reversibility of the dynamics and Lemma [23 the 
random vectors (Xo, , R^o (0, s)) and (X°, Xo, Rxo(0, s)) have the same law. Hence 

ds a.e., a.s., [/(Wo)|X°] =u(s,X°) and [i?xo(0, s)V/(Xo)|X°] = Vu{s,X°). (66) 

For s such that Equation (I 66 II holds, one deduces that 

E^„ [v/(Xo)-i?5o(0,s)aSf’A(X°)] =E^„ [E^„ [Rxo{0,s)Xf(Xo)\X°] ■ dlFxiX^)] 

= E,„ [Vu(s,X°).aSEA(X“)] 

= f Vu(s,x) ■ dxFx(x)e~^^^^dx 

= f uis,x) (XV{x) ■ dlFxix) - X ■ dlFxix)) 
ilR'i 

= E,„ [/(Xo) (VF(X“) • d°xFx(X°) - V • aSE'A(X°))] , 

where we used Lemma m below with v{.) = u{s, .) which is in H^(e ^^^^dx) for the last but one 
equality. □ 

Lemma 36. Let v be a function in H^(e~'^^^^da;). Then 

[ Xv(x) ■ d^F\{x)e~'^^^^dx = f v(x) (XV(x) ■ dxFx(x) — X ■ dxFx{x)) e~^^^^dx. 

jR'i jR'i 

Proof. Let X"(®) = x(®/’^) where X is a smooth, [0, l]-valued, cutoff function such that x(®) = 1 
for |a:| < 1 and x(*) = 0 for |a;| > 2. By integration by parts, one gets 

f Xn(x)Xv(x) ■ dxFx(x)e~^^'^^dx = f Xn(x)v{x) (XV(x) ■ dxFx(x) - X ■ dxFx(x)) e~^^'°^dx 

JRfi J'R<i 

- f v(x)Xxn(x) ■ dxFx(x)e~^^'^^dx. 

VlR'i 

The result then follows from Lebesgue’s theorem by taking the limit n —>■ 00 , using the fact that 
VF G L^(e“'^^^^dx), dxFx G L°° and V • dxFx G L^(e“'^^“^^da;) from Assumptions (Pot)-(M) 
and (Drift)-(i)-(M). □ 

By combining (1621) and (16511 . one gets: for any t > 0, 

aSE.„ If(Xt^)] = r E^„ [/(Xo) (VF • d'iFx - X ■ dlFx) (X°)] ds, 

Jo 

where, we recall, the process (X°)t>o satisfies with an initial condition Xo being distributed 
according to the equilibrium measure tto. Taking formally the limit t —^ 00 , one obtains the classical 
Green-Kubo formula discussed in the next section. 

4.2 Infinite time Green-Knbo formnla 

Theorem 37. Consider the process (X°)t>o satisfying © with an initial condition Xo being dis¬ 
tributed according to the equilibrium measure tto. Assume that Assumption (Poinc( 7 ;)) holds for 
some positive rj. Then, for any f G (e~^dx), A 1 —>■ fd-xx is differentiable at \ = 0 with 
derivative 

dl f /dTTA = r E,„ [/(Xo) (VF • d°xFx - X ■ d°xFx) (X°)] ds. (67) 

JlRrf Jo 
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Let us recall some results and notation from Section [Ml The generator of the process (X°)t>o 
is Co = —W • V + A. The generator of the process {Xt)t>o is Cx = Fx ■ V + A = Co + T\ where 
7a = {Fx + W) • V. The domain of the operator Co is 

V{Co) = £ Lo(e~'^*'“^^da;) n BI^(e~^*'“^^da;), Cov £ Lo(e~'^^“’^d2:)| . 

For any / £ there exists a unique function g = —inL§(e“'^*'^M 2 :)nEI^(e“'^^“^^da:) 

such that 

Vu £ Lo(e~'^*'^^da:) n lHI^(e~'^*'^^da:), f Vg{x) ■Vv{x)e~'^^^^Ax = f f{x)v{x)e~^^^^dx. 

Jit’i' Jit<i 

Let us start with a lemma which is a consequence of the results of Section m\ on the long-time 
behaviour of 

Lemma 38. Let us assume that Assumption (Poinc(7;)) holds for some positive p. Let us introduce 
the semigroup Ft associated to the Markovian evolution Q.' for any f £ L,o{e~^^^^dx), 

Ptfix) = m{Yn) 


where (YF)t>o satisfies (I12II . ITe then have the following Laplace inversion formula for the opera¬ 
tor To for any f in Lg(e-'^(">da:), 


nc 

-Fo^f = 

Jo 


Ptfdt. 


( 68 ) 


Proof. From Proposition!^ we know that u(t, x) = Ptf{x) = ¥,[f(Yf^)] is well defined in Lo(e 

and satisfies the partial differential equation (I23II in the sense of Definition [6] From Proposition 1111 

(since f{x)e~''^'^'>dx = 0) 


Vt 


> 0, ||7i/|lL2(e-V(x)d^) < e ’’ ll/llL 2 (e-V(x)d^) . 


(69) 


This shows that Ptf dt is well defined in L§(e '^^^'^dx). Moreover, an adaptation with A = ip = 0 
(and thus Co = C = 0) of the first energy estimate in the proof of Proposition (3 shows that 


f f |VPt/(a:)|^e < oo. 

Jo Jm.’I 


(70) 


Therefore, Ptf dt £ EI^(e ^^^^dx). 

From DefinitionjSl for any test function v £ L,'^(e~^^^^dx) D E[^(e“'^^“^^da:), 

f f{x)v{x)e~^^^'’dx= f Ptf{x)v{x)e~^^^'’dx+ f f VPsffx) ■'Vv{x)e~^^^^dxds. 

JlRa Jv’- Jo JR'* 

From Equation (I69II . limt_>oo fj^d Ptf(x)v(x)e~^^^^dx — 0 and thus, from (I7UI) . for any test function 

V £ L^(e-'^'^^^dx) nM^e-'^'^^^dx), 

f f {x)v{x)e~^^^^dx = f X ( f Ptf{x)dt^ ■ Vu(a;)e~'^^“^^dx. 

aRd J-Rd \Jo J 


This concludes the proof. 


□ 
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We recall that Do the orthogonal projection from L^(e onto Lo(e (see ®)- 

We can now give a different expression for the right-hand side of (in3- From Lemma [HHl applied 
to the constant 1, one has ■ d^Fx — V • dxFx){x)e~^^^^dx = 0. Then, using successively 

this equality, the self-adjointness of Ft in \?{e~^^^''Ax) (which is a direct consequence of m), 
Equation (1681) and Lemma [36l one has. 


[f(^o)(VV-d°xFx 

- V- 

d°xFx){X°)] df 

Jo Jtt'^ 

E[fix){VV-d°x 

Fx-^ 

^■d°xFx){Yn] e-^^^^dxdt 

Jo JRrf 

E[nof{x){\7V- 

d°xFx 

-V-d°xFx){Yn] e-'^^^^dxdt 

Jo 

nof{x)Pt{vv ■ 

dlFx- 

- V • 9A7A)(2;)e~'^*'“^^da;dt 

Jo JR'^ 

(Ptno/(a:))(vy 

■dlFx 

— V • i9A7L\)(2:)e~'^^^^da;dt 

- f (^1 

7R<i 

7ino/)(a:)(VR. 

d°xFx 

-V ■d°xFx){x)e-^‘-^^dx 


= -/ dxTx[{€.o^^of)]{x)e 

where dxT\ stands for the operator dxFx ■ V (consistently with the definition (1441) of 7a). As a 
consequence, proving Equation (1671) boils down to proving: for any / G dx), 

dll Uof{x)dnx{x) = -f aSrA[(£o^no/)l(a;)e-'"("W. (71) 

JRd, 

We are now in position to complete the proof of Theorem Ezl 

Proof of Theorem \3T\ Recall that for A small enough, the operator I + (TaTIo^^IIo)* is invertible 
from lJ^{e~^^^'^dx) to itself with bounded inverse (see Lemma fTTl) . For such a small A, one has the 
equality 

(7 -h (TaTIo ^Ho)*)-' = 7 - (TAro'^Ho)* + Tlx, (72) 

where ('by (H^ 'l the remainder 77.x = (7-|-(7a71o ^no)*)~^((7A71o ^IIo)*)^ has a norm from h^{e~^^^''dx) 
to itself of order Cl(A^). Thus, by the analytical formula for tta obtained in Lemma fTTl 

[ f{x)d'K\{x) - f f{x)d'Ko{x) = -f [{T\jCo^^o)*l]ix)f{x)e~^^'^^dx+f [7?.Al](a;)/(a:)d7ro 
jRa jRa jRa jRa 

= -/ [(rA7:o’'no)/l(a:)e-''(">da:-hO(A"). (73) 

Since, according to (Drift)-(i), jg bounded by C for A G (0, Ao], one has, by Lebesgue’s 

theorem, 

lim i [ T\Co^Tlofix)e~'^^'^^dx = [ {dxT\)jCo^^ofix)e~^^'^^dx. 

A->0 A 7]gd J]gd 

Dividing dZSl) by A and taking the limit A —>■ 0, one concludes that A i-^ f^d fdTvx is differentiable at 
A = 0 and (1711) holds. □ 

Combining the previous result with (1651) . we obtain the following corollary. 
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Corollary 39. Let f G L^(e be a Lipschitz function. Also assume that the Assump¬ 

tions (Poinc(77)) (for some positive p) and (min Spec) are satisfied. Then, one has 

d° ^ [v/(^o) • Rxo (0 . s)9S^a(X°)] ds. 

5 Long-time convergence of the estimators (|TU1) and ( fTTI) 

5.1 Statement of the main result 

Let us study the long-time behavior of the two estimators (HSl) and m- 
Theorem 40. Let / : — >■ R fee a function such that V/ is bounded. 

• Assume the existence of p > 0 such that either V is p-convex or Assumptions (Poinc(T;)) 
and (Spec(? 7 )) hold. Then A i—> i /(V^)ds is differentiable at X = 0 with derivative 
i/o‘V/(X“).T.ds and 




(74) 


• Assume either that V is p-convex for a positive constant p > 0 and E|Vo| < -boo, or that 
there exist p > 0 and p G (l,oo] such that (Poinc(r;)) holds, Xo is distributed according to 
a measure having a density in IX{e~'^^^'^Ax) with respeet to e~'^^^^Ax and (Spec(p)) holds 
for some p < p{p — l)/p (with the convention p{oo — l)/oo = p). Then VA G [0, Ao], Vfe > 0, 
E|/(X/')| < -boo, A 1 —^ E[/(X/')] is differentiable at X = 0 with derivative 9°E [/(V/')] = 
K[Vf{X°).Tt] and 

(75) 


hm dlE \f{X^)] =dl{f fAw^ . 
t^oo L J \J^d J 


Remark 41. When tto is assumed to satisfy a logarithmic Sobolev inequality with constant p, which 
is stronger than the Poincare inequality (Poinc(77)), then the second statement still holds as soon as 
Vo is distributed according to a measure having a density in U'{e~^^^'^Ax) with respect to e~^^^''Ax 
for some p > 1 and (Spec(p)) holds for some p < p, because of the hypercontractivity property of 
the semi-group associated with ensured by the Gross theorem. 


5.2 Long-time behaviour of (X°,Tt 


t t>o 


To prove Theorem 1401 one first needs to know the long-time limit of the trajectory and its tangent 
vector. We more generally consider (V®,r°)t>o solving 


dVf = -XV{X?)At + V2AWt 
AT° = {dlF>,(X°) - VV(V°)T°) dt 


(76) 


with (Xq,Tq) any initial condition independent from the Brownian motion (VKt)t>o. To write 
conveniently the long-time limit of {X^,T°), we will run time backward and use Lemma [22] about 
the link between the forward resolvent and its backward counterpart. 

Lemma 42. Under Assumptions (Poinc(p)) (for some positive p), (min Spec) and (Conv), the 
process (V°,T°)t>o converges in law as t goes to infinity to the couple 


(yo,J^ R^{0,t)dlFxiYt)At 
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where {Yt)t>o follows the dynamics ([2]) with with Yo distributed according to e Moreover, 

the law V of this couple is invariant by the dynamics and ergodic for this dynamics: for any 
test function c/p : x R'^ R in L^(V), for V-a.e. deterministic initial condition {Xq,Tq), 

lim i [ g}{X°,T°)ds = [ ip{x,T)dV{x,T) a.s.. 

t^oo t Jq JgiJxKa 


Proof. The integral Ry (0,t)dxFx{Yt)dt is almost surely well defined, from Lemma [24l and the 
bouirdedness of dxF\. To prove Lemmawe are going to use a time reversal argument. 

For to > 0, we construct a coupling of the trajectory (Xt)t>to with another process (Xt°)t>to 
following the dynamics m, but being at equilibrium. Denote by qtg the density of the dis- 

tribution of X^^ (which exists by Lemma [2|, and define ptg = —■ Let U and be 

mutually independent random variables which are independent of Xq and of the Brownian mo¬ 
tion (Wt)t>o driving (X°)t>o, such that U is uniformly distributed over [0,1] and, when qtg yf e~^, 
(tg is distributed according to — qtg(x))~^dx, C being a normalization constant does 

not need to be defined when qtg = e~^). We define the position of the process {Xt°)t>tg at 
time to by Xt^ = ^tg'^u<ptg{x° ) + Cto'^u>ptg{x° )> which is distributed according to tto. One 

has P(xt“ 7 ^ ^°o) = |lkto(a;) - ||Li(dx)- For t > to, let ixl°)t>tg evolve according to the dy- 

namics m with Brownian motion (Wt)t>o. Notice that (xt+to)t>o has the same law as the process 
at equilibrium {Yt)t>o introduced in the statement of Lemma l42l Moreover, {Xt°)t>tg is such that 

P(Vt > to,xl° = Xt) = 1 - ^Wqtgix) - e“'^^^^||iLi(dx). 

From an easy adaptation of Proposition 1201 one has on the one hand 

(X?, r°) = (^X°, R^O (0, t)T° + R^g {s, t)d°xFx{X°)ds 

On the other hand, for 0 < to < t, we have the equalities, by using successively the time transla¬ 
tion s —>■ s —to, the change of variable u = t — to — s, Lemma E^ ( using the notation, for u G [ 0 , t —to], 
Zu = Yt-tg-u) and the time reversibility of the dynamics ([ 2 |: 


(x\\ Rx^o{s,t)dlFx{x\°)ds^ = (Yt-tg,l^ Ry is, t-to)d°xFx{Ys)ds'j 

= (Yt-tg, j Rvit-to-u,t-to)dlFx{Yt-tg-u)du'j 
Zo,J^ ° Rl{d,u)dlFx{Zu)du 
Yo,J^ R^iO,s)dlFxiYs)ds 


where = stands for the equality in distribution. As a consequence, for any bounded Lipschitz 
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function 95 : x R'* —>■ R (with Lipschitz constant Lip(95)), for t > to > 0 


E 


< 


+ 


+ 


< 


(yo, JJ' R^{0,s)d°xF^{Ys)dsj - v{X?,Tt°) 


Yo,J R^{0,s)dlF>.{Ys)dsj -^{Yo,J R^{0, s)d’iF^{Y,)ds 

Rx^o{s,t)dxFx{xi°)ds^ 

Rxo (0. t)ro° + r° Rxo (s, t)d°xFx (X°)ds + {s,t)d°xFx (x^" )d, 

^ (xl°,Rxo{0,t)T° + £ R^o{s,t)dlFx{X°)ds + I' R^t, {s,t)dlFx{x\°)ds^ - ^(X°,r°) 
V>(Yo,j^ RY{0,s)d°xFx{Ys)ds'j -^(yo,J^ ° R^{0,s)d°xFx{Ys)ds 


+ E 


(2||(/j||L=o(n)) A ( Lip((p) 


Rxco{to,t) (^Rxo{0,to)TS + Rxo{s,to)d°xFx{X°)ds 


+ 2M^¥{X?^^xll)- 


(77) 


Notice that we used the semi-group property of Rxo to obtain the last but one inequality. The first 
term in the right-hand side converges to 0 as f —>■ 00 by Lebesgue’s theorem. A direct adaptation of 
Lemma [24l shows that Rxo {to, i) goes to 0 as t goes to infinity, yielding from Lebesgue’s theorem that 
the second term in the right-hand side of (Ezl) goes to 0 as t goes to infinity. The third term in the 
right-hand side of can be rewritten as 2\\>p\\^¥{X2„ ^ Xto) = ll</’l|L~(n) - Pto{x)\\i^i 

and thus goes to 0 as to goes to inhnity, by Corollary [121 Letting t —^ 00 and then to 00 in (|Z3, 
we conclude that the couple (X°,T°) converges in law to {Yo,f^ Ry{ 0, s)dxFx{Ys)ds) with law V. 

To check that V is invariant by the dynamics (Eg. we denote by {Vs)s>o the Markov semi-group 
associated with this dynamics. One has E[Ps(p(A°,T°)] = E[(p(A)Ys,TtYs)] where the right-hand 
side converges to ip{x,T)dV{x,T) as t —>■ 00 and 'Ps‘p{x,t) = E[(p(y/,TJ'’’’)] with 




T iX,T 

t 


= x- [ Vl/(yY)ds -f V2Wt, 

Jq 

= T+ f {d°xFx{Yn - vV(yY)r;’") ds. 

^0 


The continuity of a: 1 —^ {Xt)t>o for the topology of local uniform convergence on C(R+,R‘*) together 
with the continuity of V^V implies the continuity of a; 1 —^ (7?v“= (s, t))s,t>o for the topology of 
local uniform convergence on C(R+ x R+,R‘^’'‘^). With the continuity of dxFx, one deduces the 
continuity of (a;, r) i->- Tf= Ry^ (0, s)r-|- Ry^i (r, s)dxFx{Y^)dr. Hence, by Lebesgue’s theorem, 
Ps^piXjT) is continuous and bounded and E[Psi/9(X°, r°)] converges to Ps‘p(x,T)dV(x,T) as 

t-s- 00 . Therefore Ps(p(x,T)dV{x,r) = (,a(a;,r)dV(a:,r) and the probability measure 

V is invariant. 

Since V is the unique invariant probability measure for the SDE (1761) . this measure is ergodic 
(see for example |27[ Theorem 3.8 and Equation (52)]). □ 

Let us deduce from the previous results the limit of j ‘p(Xg) ■ T°ds where 95 : R'* —>■ R'* is 
measurable and bounded. 
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Lemma 43. Assume the existence of rj > 0 such that either V is rj-eonvex or Assumptions 
(Poinc(77)) and (Spec(77)) hold. Then |r|dV(a;, t) < oo (where the probability distribu¬ 
tion V has been introduced in Lemma and for any function tp : ^ measurable and 

bounded, j • Tfds eonverges a.s. to fj^d^Rd ' T"dV(a;,r) as t ^ oo whatever the choice 

of the initial condition {Xq,Tq) independent of the Brownian motion (Wt)t>o- 

Proof. Notice that if V is ?;-convex, Assumptions (min Spec), (Conv) and (Poinc(77)) hold. 
In addition, Assumption (Spec(77)) implies Assumptions (min Spec) and (Conv). Therefore, 
the conclusion of Lemma 32] holds under the two classes of hypotheses considered. Since dxFx is 
bounded, one has 

/ |r|dV(z,T)< rE\\R'^{0,t)\\dlFx{Yt)\\dt<C rE[\RY{0,t)\]dt, 

JRd^^Rd Jq L J Jq 

where the right-hand side is hnite by (I58II when V is ry-convex and by Proposition 1271 otherwise. 

In case the law of {Xq,Tq) is absolutely continuous with respect to V, the result of Lemma |43] is 
then a direct consequence of the ergodic property of the process {Xt ,Tt)t>o stated in Lemma HD 
To extend this result to more general initial conditions, we proceed as follows. By Lemma 3] 
the law of Xi is absolutely continuous with respect to tto which is the marginal law of the d first 
coordinates for the ergodic measure V. Let dVrix^xir) denote a regular conditional probability 
distribution of the d last coordinates given the d first ones under V and T° be a random vector 
independent of {Wt — VKi)t>i with conditional law given X^ equal to (''")■ Let for t > 1, 

ft° = Rxo{l,t)f? + Rxo{s,t)d‘iFx{X°)ds. Then dft° = (^d°xFx{X?) - V^V{X?)ft°) dt so that 

(X°,f°)t>i solves (Eg starting from (X°, T®) the law of which is absolutely continuous with respect 
to the measure V ergodic for this stochastic differential equation (see Lemma 3U- As a consequence 
j (p(A°) • Tsds converges a.s. to fj^^d^md ' 'rdV(a:,r) as t —>■ oo. Now, by an adaptation of 
Proposition l2UI one can check that for t > 1, Tff—Tff = (1, t){T^—T^) = Rxo (0, t)Rxo (1,0)(Ti’ — 

ri") so that 

i p{X°) ■ T°ds p{X°) ■ f°ds ip{X°) ■ T°ds 

+ p{X°) ■ Rxo(0, s)dsRxo(l, 0)(Ti° - fi°). 

The proof is completed by noticing that this quantity converges a.s. to 0 (for the second term in the 
right-hand side, this is a consequence of the boundedness of ip and of the almost sure estimate ®). 

□ 


5.3 Proof of Theorem 140! 

We are now in position to prove Theorem 1401 


Proof of Theorem \40\ Let us start by a preliminary result concerning the integrability of the random 
variable |V/(Ao)| \R^o{0,t)dxFx{Xt)\ dt. By the boundedness of V/ and dxFx 


Etto 


■ roo . . 

/ \XfiXo)\\RloiO,t)d°xFx{X°)\dt 
Jo ' ' . 



Etto [|iixo(0,t)|]dt. 


where, we recall, the subscript tto in E.,ro indicates that Xq is distributed according to ttq. If P is 
77 -convex, using (l58ll . almost surely, 

|i?xo(0,t)| < e 
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If Assumptions (Poinc(77)) and (Spec(77)) hold for some positive rj (notice that (Spec(p)) for 
P < with p £ (l,+oo] implies (Spec(77))), then, by Proposition 1271 . 

E.0 [|i?jfo(0,t)|] < 


for some positive constant C. Hence, in all cases. 




POO I 

/ \yf{Xo)\\RloiO,t)dlFxiX?) 

Jo ' 


dt 


< + 00 . 


(78) 


Let us now prove the first statement of Theorem 1401 By (I63II . the boundedness of V/ and 
Proposition 1191 Lebesgue’s theorem implies that A i—>■ j Jp* f{X^)ds is differentiable at A = 0 with 
derivative j V/(A°) • Tsds. By Lemma|43l \ Xf{X°) ■ Tsds converges a.s. to V/(a:) • 

rdV(a;,T) = [V/(Ao) • /p°°-R^o (0, s)9°i<A(A's )ds]. The proof of (iTil) is then completed by the 

following computations: 


E-n-p 


POO 

V/(Xo)-/ Rlo{0,s)d‘lFx{X°)ds 
Jo 


= ^°°E,p [v/(Xo) ■ R5 o(0,s)95Pa(A°)] ds 

=d° [ f{x)d-Kx{x) 


(79) 


where we used Fubini’s theorem and (I78II for the first equality and Corollary 139 1 for the second one. 

Let us hnally deal with the second statement of Theorem|40l By (I63II . the boundedness of V/ and 
Proposition 1191 it is enough to check that E|/(X°)| < +oo to deduce that VA € [0, Ao], K\f{Xt )| < 
+00 and A i—>■ E[/(X^)] is differentiable at A = 0 with derivative 9°® [/(A^)] = E[V/(X°) • Tt]. 
When satisfies a Poincare inequality, then according to [4] and the references therein, since 

/ is a Lipschitz function, there exists a positive e such that Jgd e~^^^^dx < +oo. Therefore, 

when the law po of Ao has a density with respect to e~^^^^dx in lJ’{e~^^^''dx), by Lemmajd] 


supE|/(A?)| < 

t>0 




dpo 


e '^C)cia; 


< + 00 . 


]LP(e-V(x)dj;) 


When V is 77 -convex, then computing \Yt\'^ by Ito’s formula, remarking that 
-2VV{Yr) • T" = -2{VV{Yr) - Vl/(0)) • Y^ - 2VP(0) • Y^ < + 


.,2 , |VP(0)|^ 




applying a localization procedure to get rid of the expectation of the stochastic integral, one obtains 


Vt > 0, E[|yt"|^] < e-''‘|a;r + 
with law po is integrable. 


-Vtl^l2 , 1-e-^* ( |VV(0)| 


+ 2d). Hence, when the initial random variable Ao 


E|At°| < [ E|W"|d^o(a:) < 

dlR<i 


dlR<i 


Y-\^]dpo{x) 


< E 


e-’J^IAoP + 


1 - / |vy(o)p 

V \ V 


+ 2 d 


< + 00 , 


and E|/(A)’)| < +00 by the Lipschitz continuity of /. 

Notice that if P is p-convex, Assumptions (min Spec), (Conv) and (Poinc(77)) hold. Moreover, 
Assumption (Spec(p)) implies Assumptions (min Spec) and (Conv). Therefore, the conclusion 
of Lemma[33holds under the two classes of hypotheses considered. The function {x, r) 1 —V f(x) ■ r 
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is continuous and the family {Vf{Xt) ■Tt)t>o is uniformly integrable by Proposition 1291 when V is rj- 

/. j, v(p-i) \ _ 

convex and since supt>Q E I jV/(Xt) • Tt I ’’p J < oo. bv Proposition 1301 in the second framework. 
Therefore the convergence in distribution in Lemma l42] yields 

POO 

V/(Xo)-/ Rlo{0,s)dlFx{X°)ds , 

Jo 

where the right-hand side is equal to 9° f(x)dTY\(x) according to (I79II . □ 


hm E (Vf(X°) ■ Tt) = E,„ 
£—>■00 


6 Numerical illustrations 

In this section, we illustrate through various numerical experiments the theoretical results obtained 
above. In Section [6T] we study numerically on a one-dimensional toy model the integrability of the 
tangent vector Tt and the sharpness of the integrability exponent obtained in Proposition 1301 In 
Section ED we illustrate the interest of the estimator m on a more realistic test case proposed 
in [34]. Finally, we investigate in Section |6]3] on the one-dimensional toy model a variance reduction 
method for the estimator (na. 


6.1 A one-dimensional toy model 


In this section, we would like to study on a simple test case the integrability of the tangent vector Tt, 
and to compare the theoretical bounds obtained in Proposition (30] with a numerical estimation of 
the integrability exponent. Let us consider the potential 

Vx € R, Vx{x) = x''‘ — ^x^ + Xx, 


where c is some fixed constant, and A G R is the parameter. For A = 0, Vb has curvature —c at the 
origin, and for c > 0, Vo is a double-well potential, with wells located at ±y^/2 and separated by 
a barrier with height (?/Vo. In particular, as c gets larger, the dynamics (|2J of {Xt)t>o becomes 
more and more metastable. 

Let us start with some explicit computation on Tt. When A > 0, the perturbative force pushes 
the system to the left. Therefore, one expects the tangent vector Tt to be negative in the mean. In 
fact, one can prove that Tt is in that case almost surely negative for f > 0. Indeed, {Tt)t>o satisfies 
the equation 

f dtTt = -l + {c-12{X/ f)Tt, 

1 ro = o, 

which can be solved explicitly, since in dimension 1, the resolvent Rxo{s,t) is given by the expo¬ 
nential Rxo{s,t) = exp (c{t — s) — 12 jJ(Xu)^du). Equation (I56II then becomes 


Tt = — exp — s) — 12 J (X°)^du^ ds < 0. 

Concerning the upper bound on the integrability exponent obtained in Proposition 1301 if the 
initial condition Xo has a bounded density, then the tangent vector is bounded in L“ uniformly in 
time for all a strictly smaller than rj/p. Here, rj is the Poincare constant of the potential Vo and p 
is the quantity 


P 


— (inf min Spec(V^ V(a;))) 


(minSpec(V^V'(a:)))^ e 
(find min Spec(V^P(a;))e“'^C)d 2 ;)^ 


(fu(12x2 


c)2g-VWd2; 

(80) 
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appearing in Assumption (Spec(77/Q)). The real number p can easily be approximated by one¬ 
dimensional numerical integration. Concerning the Poincare constant rj of Vo, let us first notice 
that the potential Vo can be written as the sum of a convex potential and a bounded perturbation 
and thus satisfies a Poincare inequality thanks to the Holley-Stroock perturbation lemma (see [T] 
Theorem 3.4.1]). The corresponding Poincare constant can be computed numerically, since it is the 
second eigenvalue of the operator L = — Vo{x)dx = e^°dx{e~'^°dx), whose first eigenvalue and 

eigenvector are 0 and the constant function 1. The numerical method we use to approximate 77 
is the following. First, notice that the spectrum of the operator L is identical to the one of L = 
g-Vb/ 2 ^gV'o /2 _ gVo/ 2 Q^^^-VoQ^^^Vo /2 ^ which is self-adjoint in the space L^(da;). The operator L 
is then discretized using a regular mesh with constant space step 5x by the infinite tridiagonal 
matrix defined by 


Mi,i = - 


— 

V 


— V((i+l/2)5x) ^ ^V(i5x) — V((i —l/2)i5x) ' 


and 




^ < 5*2 


(with Mij = 0 whenever \i—j\ > 1). We consider the restriction to a finite set of indices (Mi,j)_jv<i,j<jv, 
which is equivalent to imposing homogeneous Dirichlet boundary conditions at a: = —NSx and 
X = NSx. These artificial boundary conditions are justified (in the limit N —^ 00 ) by the fact that 
the eigenvectors of L go to 0 at infinity. Since the matrix {Mij)-N<i,j<N is a nonpositive symmetric 
matrix, one can successively compute its first eigenvalues by the inverse power method, using at 
each step a projection on the orthogonal of the eigenvector which have already been computed. We 
checked that the numerical approximation obtained for the second eigenvalue is converged when 
(5a: —>■ 0 and N ^ 00 . The graphs of numerical approximations of both p and the Poincare constant 
are plotted on Figure[T] In particular, for a curvature constant c located left to the intersection of the 
two curves (approximately c < 0.86), Proposition 1301 ensures that Tt is bounded in L^, uniformly in 
time. Also, for curvature constants such that p is less than half the Poincare constant (correspond¬ 
ing approximately to c < 0.50), Tt is bounded in L^, and thus has a bounded variance uniformly in 
time. On Figure [2] we plot the critical exponent rj/p such that, according to Proposition 1301 Tt is 
in L“ for a < pip. 

Let us now explain how we estimate numerically the integrability exponent a such that Tt actually 
is in L“. This is done by computing the tail of the empirical cumulative distribution function of Tt- 
We simulate 10® independent realizations of the process {Tt,Xt), starting from ^/c/2 (that is, at the 
bottom of the right well), up to the time t = 40, at which the systems seems to be at equilibrium. 

On Figure [S] we plot in logarithmic scale the tail of the empirical cumulative distribution function 
of those N = 10® independent realizations (Tt )i=i,...,jv, namely 


1 1 

[0,oo) 9 0 :^ 


N ■ 


with curvature c being respectively 2, 3, 4 and 5, from bottom to top. Linear regression in those 
four cases gives the following slopes: 


c 

2 

3 

4 

5 

slope 

-3.09 

-1.95 

-1.29 

-1.12 


We have checked that the results are the same for f = 40 and for t = 80. Note that an integrable 
random variable corresponds roughly to a slope less than —1, and a square integrable variable cor¬ 
responds to a slope less than —2. We also plot on Figure [2]the empirical integrability exponent for 
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Figure 1: Poincare constant of the measure e (dashed line), and parameter p defined in H80I1 (solid line), as 
a function of the curvature c. 



oi-^^^^^-1 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Curvature at the origin 


Figure 2: Blue: theoretical lower bound of integrability for Tt, according to Proposition [30] as a function of the 
curvature c (In fact, Proposition [3D] yields boundedness in L“ only for a > 1, corresponding to c < 0.86) ; Red : 
Numerical estimation of the integrability exponent a (for t = 40). 
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Figure 3: Logarithmic plot of the empirical cumulative distribution function of 10® independent realizations 
of Tt (for t — 40), with parameter (from bottom to top) c = 2, c = 3, c = 4, c = 5. 


different curvatures between 0 and 3. We observe that the results are in accordance with Proposi¬ 
tion ESI the theoretical lower bound is indeed smaller than the effective integrability exponent. 

For a curvature larger than 3, the tangent vector Tt at time t = 40 does not seem to be of finite 
variance. This raises the question of appropriate variance reduction technique to be used in order to 
use the estimators (nsi) or m- We will investigate in Section [6.31 a hrst idea that could be used in 
this one-dimensional situation. Further studies related to this problem will be the subject of future 
works. 


6.2 A many particle system 

In this section, we consider a more complex problem introduced in El]j and motivated by exper¬ 
imental studies of colloidal particles in optical traps. Let us consider Xt = (F)^,- • ■ ,Yt^) where 
(F/)i=i.., ,N € (R^)^ are the positions of N two-dimensional particles evolving according to 

N 

dv; = -KYt’-dt + VC/(y/ - y/)dt -b Aeidt -b dwi, l<i<N (81) 

with K > 0, A € R, ei the normed vector directed along the first coordinate, and U{x) = /\x\. 

The particles undergo a quadratic confining potential near the origin with strength k, a repulsive 
interaction given by U, a shear in the x-direction with strength A and a thermal noise. 

We study the case of Af = 10 particles with repulsion range F = 25 and attraction intensity k = 10 
corresponding to the parameters studied in |34| . For those parameters, at equilibrium, particles are 
gathered around the origin. At A = 0, no particular direction appears in the dynamics, and the 
equilibrium measure is invariant with respect to rotations around the origin. 

One wants to study the effect of shearing on the symmetry of the invariant measure. This 
symmetry can be measured by the empirical covariance $ of the particle system, dehned by 

i=l 
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Figure 4: Sensitivity of the covariance with respect to the shear. 


where Yt = and for k G {1, 2}, i Y-^’’^. One is interested in computing the 

derivative 8° /(b, 2 )N 4>d7rA. 

On Figure|4l we plot the confidence interval obtained for the expectation E[9A4’(-Yt^)], with N = 
10® independent simulations, as a function of the time t. The dynamics m has been simulated 
using an explicit Euler-Maruyama scheme with time step St = 10“®, and the expectation has been 
calculated through the Monte Carlo approximation 

1 

E[9S4>(Xt^)] = E[Tt • V<I>(X“)] V<1>(X?-*) 

i=l 

where the (X^’' are independent simulations of the Euler-Maruyama discretization of the 

dynamics 0 ruling the evolution of {Xt,Tt). As in [34], we observe that the correlation function 
increases as a function of time, before reaching a plateau. We have checked that similar 
results are obtained using a finite differenciation instead of the simulation of the couple (X°,Tt). 

6.3 Particle merging 

As mentioned in Section [6. II in some situations, the variance of the tangent vector may become very 
large (or even infinite) which means that the estimators (IIOII and (Hill become ineffective. Therefore, 
it is desirable to introduce variance reduction mechanisms. We explore in this section a first idea 
in the simple one-dimensional test case of Section 16.11 Extensions and further variance reduction 
techniques will be the subject of forthcoming works. 

A first simple idea to reduce the variance is to replace the tangent vector Ta in the estimator of 
E[V/(X°) • Taj by its conditional expectation given A®. This corresponds in practice to replacing 
the tangent vector of particles which are at the same position at a given time s by the average of 
their tangent vectors. Then, the particles evolve again following the dynamics 0. We refer to this 
procedure as “particle merging”. In practice, with this naive procedure the probability to observe 
two particles at the same position is zero, in dimension larger than one. A first simple practical 
way to implement this technique is to introduce small subsets of the configuration space, and to 
merge particles which are in the same subset, which of course reduces the variance but introduce 
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some bias. The merging can be performed in a much efficient way and in larger dimensions, by 
correlating the particles, see e.g. m- This will be the scope of future work. Before studying the 
interest of particle merging on the simple case of Section 16.11 let us first state the theoretical result 
which justifies the use of this approach. 

Lemma 44. Assume (min Spec). For s > 0, let {Tt)t>B be the solution to 

[ ^ = dlF^{X°) - X\X°)ft , for all t > s, 

I f. =E[T,|X°]. 


Then Vt > s, ft = E[Tt|X°, (Wr — Ws)rg[s.t]]. Assume moreover that / : R'^ — >■ R'^ is a Lipschitz 
funetion belonging to U’{e~^^^''dx) for some p G [1, oo] and that the initial eondition Xq to (fT|l 
admits a density with respect to ttq belonging to Lp-i (where, by convention, = oo 

if p =1). Then, for each t > 0, f(Xt) is integrable for A € [0, Ao], A i—E[/(Xt^)] is differentiable 
at X — 0 and 

(E[f{Xt^)]) = K[Vf{X?) ■ Tt] = nvf{x?) ■ ft], for each t > s. 

This Lemma shows that if, at a given time s, the particles at position X^ replace their current 
tangent vectors by an average of these tangent vectors, and then follow the dynamics (O for t > s, 
the estimator m is still consistent. 


Proof. By Lemma [23l Assumption (min Spec) ensures that Tt is integrable for each t > 0. In view 
of the equality (I56II and using the semigroup property (I55II of Rxo, one gets that for t > s >0, 


Tt = Rxo{s,t)Ts+ [ Rxo(r,t)dlFx{X°)dr. 

J S 


Since (Al°)rg[s and therefore {Rx°{'ef))re[s,t] are measurable with respect to the sigma-field gen¬ 
erated by Ag and (Wr — Ws)rg[s,t], one deduces that 


E[Tt|A°, (Wr 


lT.),g[,.q] = i?;,o(s,t)E[r,|A°, (Wr 


^s)rg[s,f]] 



Rxo{r,t)dlFxiX°)dr. 


The independence of {Xf,Ts) and [Wr — Ws)rg[s,t] implies that E[rs|A°, (VIA — Ws)rg[s,t]] = 
E[Ts|A°]. Since, by an adaptation of Proposition 1201 

ft = i?jfo(s,t)E[T.|A°] -t r Rxo{r,t)d°xFx{X°)dr, 


one concludes that ft = E[Tt|A°, (Wr — l^s)rg[s,t]]- 

If the initial condition Ao to o admits a density with respect to tto belonging to Lp-i (e 
then so does A° for each t > 0 by Lemma [d] When / : R'^ —>■ R'* is a Lipschitz function belonging 
to U‘{e~'^''^'^dx), the integrability of E[/(Ai^)], the differentiability of A i—>■ E[/(At^)] at A = 0 and 
the equality 5° (E[/(A^)]) = E[V/(A°) • Tt] are deduced from an adaptation of the beginning of 
the proof of Theorem m Now, for t > s, 

E[V/(A°) • Tt] = E[E[V/(A°) • Tt]X°, (Wr - lP.).g[.,t]]] 

= E[V/(A°) • E[Tt]X°, (W. - W«),g[,,tj]] = E[V/(A°) • ft]. 


□ 
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Figure 5: Empirical average and 95%-confidence interval for the Monte Carlo estimator of 9AE[/(-^t^)] as a 
function of time, / being a smooth approximation of the the indicator function of R+: f{x) = ^ + -^ arctan(10a:). 
The estimator is built with 10® realizations. Green: with merging (10® independent relizations of 10® interacting 
particles); Red: without merging (10® independent relizations). 


To test the interest of this approach, we consider again the setting of Section 16.11 with c = 
2.9 (which corresponds to case where the variance of tangent vector Tt, at t — 40, is very large, 
see Figure [2). The merging procedure is done as follows: a uniform mesh with step size 0.04 is 
introduced, and, every ten timesteps, the tangent vectors of particles which are in the same bin 
are replaced by an average of these tangent vectors. On Figure [S] we observe that this procedure 
divides approximately the variance by 4, while introducing a bias which is sufficiently small so 
that the confidence interval of the simulation with merging is included in the confidence interval of 
the simulation without merging. Figure [2 then gives more quantitative estimates of the variances 
of these two simulations (with and without merging), as a function of time. We have observed 
numerically that large values of Tt become very unlikely with the merging procedure: using 10® 
independant realizations of 10® interacting particles over the time interval (0,10), we did not observe 
any realization of Tt with absolute value larger than 3 (compare with what is reported on Figure [JJ. 

A Alternative bounds on the density of 

In this section, we would like to present a few results which can be obtained under the assumption 
Assumption (V). The function V is of class and satisfies 

Cv = sup (2AF(z) — IVF(a:)!^) < +oo. 

Note that simple assumptions on the quantity av{x) = 2AF(a;) — |VF(a;)|^ can give strong 
results on the equilibrium measure e~^^^^dx. For instance, if av{x) goes to —oo at inhnity, then 
the equilibrium measures satisfies a Poincare inequality (see for example the appendix in |32|1. 
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Figure 6: Empirical variance for the Monte Carlo estimator of <?AE[/(-^t^)] as a function of time, / being the 
indicator function of 1R+. Green: with merging ; Red: without merging. 


A.l Bounds on the density of 

Proposition 45. Consider the setting and the notation of Lemma^and let Assumption (V) hold. 
Assume that the measure e^^dfio can be written as 

= f{x)dx + dv, 

where f is some function in L*’(da;) with p G [1, 2] and v is some finite measure on Then, for 
any t > 0, pit is absolutely continuous with respect to e“2^G)dx with 


dfit 


e-i^Odx 


< Ce'- 


L2(da:) 


a 


+ 


^(l/p-l/2)d/2 ' ^d/4 


(82) 


Proof. Let i/i : —>■ R be a bounded measurable function and recall the formula 

E[^{X?)] = E [r/>(Xo + /„‘(2Av-|vv|=)(Xo+V2U/.)d.J 

obtained by the Girsanov theorem, see Equation m- 

If C is an upper bound for i(2AF — |VF|^), and if one assumes if > 0, one obtains 


J ifdpt < e'^*E [V'(^o + 


e 

(47r; 


Z' ip{y)e { f e^^^'°'>e d/ro(x)^ di/. 

JRd yjRd / 


SO that pit ^ e i'^^^^dx with a Radon-Nikodym derivative —— satisfying 


dpt 


e-2^C)dx 


< 


2 * Tt = {f * jt + n * ■yt) , 


48 






























1 

where * stands for the convolution product, and 7 t(a:) = ^^ 77^^776 denotes the centered Gaus¬ 
sian density with covariance matrix 2tld- One concludes that (I82II holds by: 

• the Young inequality ||/*7t||L2(d:i:) < ll/l|LP(d 2 ,)|| 7 t|k<i(d:i:) where 1/p + l/q = 3/2 {p,q G [l,oo]) 
and the heat kernel estimate || 7 t||LQ(dcc) ^ Cqt \ <iJ ^; 

• the estimate |k * 7*11^2(d^) < Ik * 7t||Li(dx) Ik * 7t||L°°(dx) < 

□ 


A.2 An additional result 


Assumption (V) can also be useful to prove the second point in Assumption (Pot)-(M) on the 
potential V. 

Lemma 46. Under Assumption (V), the function VY is in 



< Cv 


Proof. Let x>i(®) = x(*/^) where X is a smooth, [0, l]-valued, cutoff function such that x(®) = 1 
for kl < 1 and x(®) = 0 for |x| > 2 . 


I • Ve 


-V{x) 


f \VV\^(x)xn{x)e = - / (xn(x)VY(x)) ■ 

JRii JRii 

= f Vxn(x) • Vy(x)e~'^^^kx + f Xn{x)AA/{x)e~'^^^'’Ax 
< f Axn(x)e“'^^'^kx + i f Xn{x)\Wf{x)e~^^'°'’Ax+]-Cv- 

Ja'i 2 I 


As a consequence. 


f |VYk(x)xn(x)e < 2 / Axn(x)e ^^"^^dx + Cv 

jR'i jR'i 


and the result follows from taking n —>■ 00 by Fatou’s lemma for the left-hand side and Lebesgue’s 
theorem for the right-hand side. □ 


B About the Assumption (Conv) 

In this section, we show that Assumption (Conv) is a natural one, since it appears as a sufficient 
condition in another related problem. 

We recall that {Yf^)t>o is defined in (US as the solution to 

Vt>0, Yt" =x - [ VV{Y^)ds + y/2Wt. (83) 

Jo 

Since {RY^{0,t))t>o is the differential of the trajectory (Yk)t>o with respect to x, we expect that 
a condition yielding long-time decay for will imply that trajectories with same noise and close 
initial conditions will eventually converge toward each other. More precisely, we are interested in 
the joint long-time behavior of the so-called duplicated dynamics {Yf^,Yf^)t>o, where x and y are 
two different initial conditions. Note here that the two processes {Yjf)t>o and (y/)t>o are driven 
by the same Brownian motion. 
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In [^, the same problem is considered for a diffusion whose diffusion matrix may not be constant. 
In that case, an example is provided, where the process does not converge to 0. 

A similar problem is considered in [3: the process is a Brownian motion reflected on the boundary 
of a domain Q,. Such a dynamics can be formally seen as a singular case of the problem we consider, 
with V = oo X Ifjc. Equation (I83II then has to be written with a local time on the boundary in 
place of W. In that case, the difference Yf^ — Y^ will converge to 0 if the domain is smooth 
enough and has at most one hole. However, it is conjectured that the same result holds for much 
more general domains. 

We will use the fact that V is such that the dynamics (1831) is ergodic with respect to the invariant 
measure 


B.l The one-dimensional case 

In the one-dimensional case, this question is especially simple, because of the order structure on the 
state space. In particular (see [23]), it can be checked that if for any a; G R, Y-i^ converge weakly to tto 
as t 00 , then the only invariant distribution of the duplicated dynamics is the image of 
by a; I— (x,x). Actually, under additional assumption, one can show that Y/^ — Y/ converges in 
mean to 0 in the long-time limit. 

Proposition 47. Assume that the dimension is d=l. If for any a; € R the time marginals of the 
process {Yjf)t>o converge weakly to no as t ^ oo and the random variables {Yjf)t>o are uniformly 
integrahle, then, for any x,t/ £ R, the process (Yff — Y^)t>o converges to 0 in L^(r2). 

According to Corollary 1121 the long-time convergence of the marginals holds for instance if the 
potential V satisfies a Poincare inequality (see Assumption (Poinc(r 7 ))). 

Proof. First, from the uniform integrability of {Xf)t>o and the weak convergence of the time 
marginals, both and converge to J^xe~^^^^dx as t oo. 

Now assume, without loss of generality that x < y. Then, from a comparison theorem, Yf^ < Yf' 
holds for all positive times, and one obtains 

E[|yr - Yt^\] = ^Yt - Yf] = ¥.\Yr] - E[y/] ->■ 0. 


□ 


B.2 A general criterion 

Proposition 48. The following facts hold true: 

1. Assume that 

yx,y G R^ (x-y)- (VH(x) - VP(y)) > - y\^ (84) 

with u : R'* —>■ R such that max(0, —ii(x))e“'^^“^^dx < oo and v(x)e~^''^^dx > 0. Then 
for all x,y £ R"^, \Yi — Y^\ converges a.s. to 0, exponentially fast at any rate between 0 
and v{x)e~^^^''dx as t oo. 

2. The exponential convergence to 0 still holds if V is convex and there exist xo G R*^ and e > 0 
such that the ineguality inf 2 ,gB( 2 ,u^j) minSpec(V^H(x)) > 0 holds. 

3. If V is convex, then the only invariant measure of the duplicated dynamics is the image 
of e~^^^^dx by X {x,x). 

Let us start with a few remarks: 
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Remark 49. 

^(1*1"- 


• The first point can be applied to the so-called Mexican hat potential V{x) = 
7 |a;p), with /3 > 0 and 7 > 0, in dimension d> 2. For this potential, one has 

{x-y)- {VV{x) - W{y)) = 2/3\x - y\^ {\xf + \y\^ - 7 ) + 2^(l®l^ - 


for v{x) = /3(4|a;p — 27 ). 

Jfld e/3(T|a:p-|3:|4)fl3. 


In addition, one has v{x)e > 0 since 

“ ")dr 

7 , 1{^=2} ^ (d-2)/+°°ri-V(^’-’-^)dr 

2 4^/+“ri-ie/3(-y-—")dr 


• Letting y ^ x in (I84II . one obtains that \/x € R'*, v(x) < min Spec(V^y(a:)). When x 1 —^ 
min Spec(V^\^(a;)) is concave, 

{x-y)- {VV{x) - VV{y)) = f\x - y) ■ VV(fe + (1 - e)y){x - y)dd 

Jo 

'>\x — y\^ f minSpec(V^V(0a; + (1 — d)y))dd 
Jo 

> \x — y\‘^ f dminSpec(V^17(a:)) + (1 — d)minSpec(V^F(y))dd 

Jo 

> i (minSpec(V^y(a;)) +minSpec(V^F(i/))) \x — y\^ 

and one may choose v{x) = minSpec(V^V^(a:)) in (1 84 II . 

• When V = V + V with V such that x 1 —>■ minSpec(V^t^(a:)) is concave and V such that x 1 —^ 
VV{x) is Lipschitz with constant S and constant outside some Borel subset A ofMJ^, then one 
may choose v{x) = minSpec(V^t^(a;)) — 2S1a{x) in (1 84 II . 


Proof. 1. One has 

d|y^- - y/|2 = -2{Yt" - Yf) ■ {VV{Yt") - vy(y/))dt (85) 

<-(o(yt") + «(y,"))|yf-y,''|^dt, 


under (IMll . Hence 

\Yt^ - <\x- y| 2 g-/o‘o(no+Hn“))d.^ 

Since, by (ISTIl . i /o(n(y/) + ?;(yi'))ds converges a.s. to 2 J^,iv{x)e~^^^'^dx > 0, one easily 
deduces the hrst assertion. 

2. When V is convex, then t 1 —^ {Yf' — YJ^\ is nonincreasing by (|8^. Now, for 2 G 5 ( 0 : 0 , f) 
and w G R'*, one has 

{z-w)-{\/V{z)-\7V{w)) > \z-w\ inf minSpec(V^y(-)) + \z - e)(u;)') . 

B(xQ,e) \Z / 
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As a consequence, 


{Yr-Yy)iVV{Yr)-VV{Yy)) > inf minSpec(VV(.)) 

B(xo,e) 


Al]\Yr-Y,^\\ 


,2|a; - y\ 

One concludes by arguments similar to the ones used for the first assertion. 

3. Let V be convex and differentiable and let a; 7 ^ y be such that {x — y) ■ (VV^( 2 ;) — W(y)) = 0. 

c+y \ _ V (3:)+' 

2 1 2 


Then V is affine on the segment [x,y] and = LlMiAlM. Por 2 e R'* \ {0} and e G R, 


V{x) + V{y) 


= V 


' 3 : + y\ ^ V(x + ez) + V(y- ez) 


(^) 


^ + f • {^V{x) - VV{y)) + o{e) 


as |e| —>• 0. As a consequence 2 ■ (V’L(x) — VV'(j/)) = 0 and W{x) = W{y). 

Let {Xt)t>o and {Yt)t>o be two solutions to the stochastic differential equation ([2]), such 
that (AojLo) is distributed according to some invariant probability measure of the dupli¬ 
cated dynamics. Since \Xt — Yt\'^ is a.s. non-increasing with t and constant in distribution, 
a.s. 1 1 —>■ \Xt — Tip is constant and therefore dt-a.e. {Xt — Yt) ■ (XV{Xt) — XV{Yt)) — 0 which 
implies W{Xt) = W{Yt). One deduces that a.s., 1 1 -^ Xt — Yt is constant. 

Now, since x 1 -^ is integrable, then V cannot be affine in some direction and for any 2 G 

R'* \ {0}, X z- {XV{x) — W(x — z)) is not constant equal to zero. By continuity of XV, one 
deduces the existence of y G R"^ and e > 0 such that Vx G B{y,e), z- {XV{x) — XV{x — z)) > 0. 
With the ergodicity of (Xt)t>o and the fact that dt-a.e. (Xo — Yq) ■ (XV(Xt) — XV(Xt — Xq + 
To)) = 0, one concludes that a.s. Xo = Yq. 

□ 


Acknowledgements 

This work is supported by the European Research Council under the European Union’s Seventh 
Framework Programme (FP/2007-2013) / ERC Grant Agreement number 614492 and by the French 
National Research Agency under the grant ANR-12-BS01-0019 (STAB). The authors would like to 
thank fruitful discussions with G. Stoltz on nonequilibrium methods and Green-Kubo formulae. 


References 

[ 1 ] C. Ane, S. Blachere, D. Ghafai, P. Fougeres, 1. Gentil, F. Malrieu, C. Roberto, and G. Scheffer. 
Sur les inegalites de Sobolev logarithmiques. Societe Mathematique de France, 2000. In French. 

[2] R. Assaraf, M. Caffarel, and A.G. Kollias. Chaotic versus nonchaotic stochastic dynamics in 
monte carlo simulations: A route for accurate energy differences in n-body systems. Phys. Rev. 
Lett, 106:150601, 2011. 

[3] V. Bally, M.P. Bavouzet, and M. Messaoud. Integration by parts formula for locally smooth 
laws and applications to sensitivity computations. Ann. Appl. Prob., 17:33-66, 2007. 

[4] S. Bobkov and M. Ledoux. Poincare’s inequalities and Talagrand’s concentration phenomenon 
for the exponential distribution. Probab. Theory Related Fields, 107(3):383-400, 1997. 

[5] K. Burdzy, Z.-Q. Chen, and P. Jones. Synchronous couplings of reflected brownian motions in 
smooth domains. Illinois Journal of Mathematics, 50(l-4):189-268, 2006. 

[ 6 ] D. Chandler. Introduction to Modem Statistical Mechanics. Oxford University Press, 1987. 


52 








[7] J. Chen, G. Zhang, and B. Li. How to improve the accuracy of equilibrium molecular dynamics 
for computation of thermal conductivity? Physics Letters A, 374(23):2392 - 2396, 2010. 

[8] G. Ciccotti and G. Jacucci. Direct computation of dynamical response by molecular-dynamics 
- mobility of a charged Lennard-Jones particle. Physical Review Letters, 35(12):789-792, 1975. 

[9] G. Ciccotti, R. Kapral, and A. Sergi. Non-equilibrium molecular dynamics. In S. Yip, editor, 
Handbook of Materials Modeling, pages 745-761. Springer Netherlands, 2005. 

[10] M. Cwikel. Weak type estimates for singular values and the number of bound states of 
Schrodinger operators. Ann. of Math. (2), 106(1):93-100, 1977. 

[11] F. den Hollander. Large deviations, volume 14 of Fields Institute Monographs. American 
Mathematical Society, Providence, RI, 2000. 

[12] M. Duflo. Random iterative models. Springer, 1997. 

[13] K.D. Elworthy and X.M. Li. Formulae for the derivatives of heat semigroups. J. Func. Anal., 
125:252-286, 1994. 

[14] D.J. Evans and G. Morriss. Statistical mechanics of nonequilibrium liquids. Gambridge Uni¬ 
versity Press, 2008. 

[15] J. Fontbona and B. Jourdain. On the long time behaviour of stochastic vortices systems. 
Markov Process. Related Fields, 20(4):675-704, 2014. 

[16] I.I. GThman and A.V. Skorohod. Stochastic differential equations. Springer-Verlag, New York- 
Heidelberg, 1972. Translated from the Russian by Kenneth Wickwire, Ergebnisse der Mathe- 
matik und ihrer Grenzgebiete, Band 72. 

[17] D. Gilbarg and N.S. Trudinger. Elliptic partial differential equations of second order. Classics 
in Mathematics. Springer-Verlag, Berlin, 2001. 

[18] M. Hairer and A.J. Majda. A simple framework to justify linear response theory. Nonlinearity, 
23(4):909-922, 2010. 

[19] Bismut J.M. Large deviations and the Malliavin calculus. Birkhauser, 1984. 

[20] M. H. Kalos and F. Pederiva. Exact monte carlo method for continuum fermion systems. 
Physical review letters, 85 (17): 3547, 2000. 

[21] 1. Karatzas and S.E. Shreve. Brownian motion and stochastic calculus. Springer-Verlag, 1988. 

[22] N. V. Krylov. Lectures on elliptic and parabolic equations in Holder spaces, volume 12 of 
Graduate Studies in Mathematics. American Mathematical Society, Providence, RI, 1996. 

[23] V. Lemaire, G. Pages, and F. Panloup. Invariant measure of duplicated diffusions and appli¬ 
cation to Richardson-Romberg extrapolation, 2013. http://arxiv.org/abs/1302.1651. 

[24] E. Lieb. Bounds on the eigenvalues of the Laplace and Schroedinger operators. Bull. Amer. 
Math. Soc., 82(5):751-753, 1976. 

[25] P. Malliavin. Integration and probability, volume 157 of Graduate Texts in Mathematics. 
Springer-Verlag, New York, 1995. 

[26] A. Quarteroni and A. Valli. Numerical Approximation of Partial Differential Equations. 
Springer, 1997. 

[27] L. Rey-Bellet. Ergodic properties of Markov processes. In Open quantum systems. II, volume 
1881 of Lecture Notes in Math., pages 1-39. Springer, Berlin, 2006. 

[28] G. V. Rozenbljum. Distribution of the discrete spectrum of singular differential operators. Dokl. 
Akad. Nauk SSSR, 202:1012-1015, 1972. 


53 


[29] T.H. Rydberg. A note on the existence of unique equivalent martingale measures in a markovian 
setting. Finance and Stochastics, l(3):251-257, 1997. 

[30] J. Tailleur and J. Kurchan. Probing rare physical trajectories with Lyapunov weighted dynam¬ 
ics. Nature Physics, 3:203-207, 2007. 

[31] R. Temam. Navier-Stokes equations, volume 2 of Studies in Mathematics and its Applications. 
North-Holland Publishing Co., Amsterdam-New York, 1979. 

[32] C. Villani. Hypocoercivity. Memoirs Amer. Math. Soc., 202, 2009. 

[33] W. von Wahl. The equation u -I- A{t)u = / in a Hilbert space and L^-estimates for parabolic 
equations. J. London Math. Soc. (2), 25(3):483-497, 1982. 

[34] P.B. Warren and R.J. Allen. Malliavin weight sampling for computing sensitivity coefficients 
in brownian dynamics simulations. Physical Review Letters, 109(25):250601, 2012. 


54 


